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DEEP HST IMAGING IN NGC 6397: STELLAR DYNAMICS 



(N 

o 

(N 
•*-> 
O 

O 

6 

in 



> 
oo 

o 
d 

(N 

> 

X 



J. S. Heyl 1 , H. Richer 1 , J. Anderson 2 , G. Fahlman 3 , A. Dotter 2 , J. Hurley 4 , J. Kalirai 2 , R. M. Rich 5 , M. 

Shara 6 , P. Stetson 3 , K. A. Woodley 1 , and D. Zurek 6 

1 Department of Physics and Astronomy, University of British Columbia, Vancouver BC V6T 1Z1 Canada 
2 Space Telescope Science Institute, Baltimore MD 21218 
3 Herzberg Institute for Astrophysics, National Research Council, Victoria BC Canada 
4 Centre for Astrophysics & Supercomputing, Swinburne University of Technology, Hawthorn VIC 3122, Australia 
5 Division of Astronomy, University of California, Los Angeles CA 90095-1562 and 
6 American Museum of Natural History, New York NY, 10024-5192 
Draft version October 5, 2012 

ABSTRACT 

Multi-epoch observations with ACS on HST provide a unique and comprehensive probe of stel- 
lar dynamics within NGC 6397. We are able to confront analytic models of the globular cluster 
with the observed stellar proper motions. The measured proper motions probe well along the 
main sequence from 0.8 to below 0.1 Mq as well as white dwarfs younger than one gigayear. The 
observed field lies just beyond the half-light radius where standard models of globular cluster 
dynamics (e.g. based on a lowered Maxwellian phase-space distribution) make very robust predic- 
tions for the stellar proper motions as a function of mass. The observed proper motions show no 
evidence for anisotropy in the velocity distribution; furthermore, the observations agree in detail 
with a straightforward model of the stellar distribution function. We do not find any evidence 
that the young white dwarfs have received a natal kick in contradiction with earlier results. Using 
the observed proper motions of the main-sequence stars, we obtain a kinematic estimate of the 
distance to NGC 6397 of 2.21^7 kpc and a mass of the cluster of 1.1 ± 0.1 x 1O 5 M at the pho- 
tometric distance of 2.53 kpc. One of the main-sequence stars appears to travel on a trajectory 
that will escape the cluster, yielding an estimate of the evaporation timescale, over which the 
number of stars in the cluster decreases by a factor of e, of about 3 Gyr. The proper motions 
of the youngest white dwarfs appear to resemble those of the most massive main-sequence stars, 
providing the first direct constraint on the relaxation time of the stars in a globular cluster of 
greater than or about 0.7 Gyr. 

Subject headings: globular clusters: individual (NGC 6397) — celestial mechanics, stellar dynamics 
— astrometry 



1. INTRODUCTION 

Globular clusters provide unique tests of star for- 
mation, stellar evolution, Galactic structure and ce- 
lestial mechanics. As inhabitants of the Galactic halo 
they probe the structure of the outer galaxy. As the 
oldest groups of stars in the Galaxy, they give a win- 
dow into the formation of the Galaxy and the early 
Universe in general — they are the foremost dig for 
Galactic archaeologists. NGC 6397 is the second clos- 
est globular cluster (after Messier 4) to the Earth, so 
many studies have focused on observations and models 
of this cluster. Globular clusters are bound groups of 
up to several million stars that all form more or less 
at the same time from the same material. A few clus- 
ters show e vidence for multip le ge nerations of st ars. In 



marks of multiple generations of stellar populations 
in NGC 6397 are much less dramatic than in othe r 
clusters such as NGC 2808 dD'Antona et al.||2005|), 

The 



flD 'Antona et al 
NGC 1851 (|Milone et al. 2008} among o thers 1. 



Lind et al.| \20ll\ and |Milone et al.| ( |2011[ ) 

re stars exhibit a bimodal abundance dis- 
In about 75% of the stars sodium, nitro- 



NGC 6397 
find that t 
tribution. 

gen, aluminum and helium are enhanced while car- 
bon, oxygen and magnesium are depleted relative to 
the remainder w hose abundances ar e simil ar to field 
stars. However, di Criscienzo et al. ( |2010[ ) find that 
the width of the main sequence allows only a small 
variation in helium or CNO abundances. The hall- 



particular abundance differenc es found by Lind et al. 
(|2011| and |Milone et al.| pOlTj ) do not significantly af- 
fect the stellar fluxes in the bands examined in this 
paper (F814W and F606W), so we cannot search for 
dynamical differences between the two populations. In 
any case given the short relaxation time of NGC 6397 
any such differences should have been washed out long 
ago. The width of main-sequence track of NGC 6397 
is not significantly larger than the observational uncer- 
tainties, supporting the approximate coevality of the 
star s and revealing a sm all binary fraction in the clus- 
ter ( |Davis et al."||2008a| ). For these reasons NGC 6397 
is an exce llent laboratory to understand stellar e volu- 
tion (e.g. [Kalirai et aLl|2008| | Richer et a L| |2008[) and 
deriv e limits on its time of formation (e.g. Hansen et al. 



2007) 



The focus of this paper is the dynamics of NGC 6397. 
NGC 6397 is the nearest cluster with a central cusp 
in its surface brightness profile, a theoretical signa - 
ture of a core collapsed cluster (Trager et al. 1995); 
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on the other hand, Heggie & Giersz (2008) argue that 
in spite of its prominent core the possibly closer clus- 
ter Messier 4 has experienced core collapse and is now 
sustained by binary burning. Core collapse is a dy- 
namical process nearly unique to globular clusters in 
which the timescale to approach thermodynamic equi- 
librium (the relaxation time) is much smaller than the 
age of the cluster and also much smaller than the evap- 
oration time (so that the cluster still exists for us to 
observe). In a core-collapsed cluster, the stars have 
achieved equipartition of their kinetic energies; mas- 
sive stars typically move more slowly than less massive 
ones. When these stars with different velocity disper- 
sions (a) occupy the same p oten tial, a signature known 



as mass segregation (see § 4.2 1 develops where more 



massive stars lie typically closer to the center of the 
cluster. Furthermore, core collapse is predicted to pro- 
ceed after the bulk of binary burning in the central re- 
gions is complete, so as the central regions of the clus- 
ter continue to lose energy to the outer regions, the 
energy transfer drives the system further from equi- 
librium, causing the core to collapse in principle to a 
singularity, but the formation and tightening of a few 
hard bina ries in the cen ter is expected to halt the col- 
lapse (e.g Spit zer|[T987 ). 

NGC 6397 has been the focus of many recent the- 
oretical investigations includin g iV — body simulations 
such a s Hurley et al. (2008) and Heggie k Giersz 
(|2009[), M onte Carlo approaches such as Giersz k Heg 
Igiel (| 2009h and le s s rece ntly the Fokker-Planck treat- 
ment of Drukier (1995). The focus of this paper is 



the dynamical observations of the cluster, so the mod- 
elling performed is less detailed than these recent in- 
ve stigations and is more in the spirit of that outlined 



Meylan k Mayor] (M9911 MM91) for NGC 6397 and 



most recently in |McLaug hlin et al.| ( |2006[ MAM06) 
for 47 Tucanae (47 Tuc). We d id compare our data 



to the recent Hurley et al. (20081 ./V— body model but 
found that the number ot the stars in the outskirts of 
the model were insufficient for a detailed comparison 
with the data. The total number of stars at the end 
of this model was about 38,000 of which about 10,000 
are useful for comparison with stars in our field — this 
resulted in the awkward situation where the statistical 
errors on the model are comparable to that of the data 
and made detailed comparision difficult. We estimate 
that a factor of three more stars would be sufficient. 
Our strategy is to fit a particular mass component of 
a multi-mas s lowered isothermal distrib ution function 
( King 1966[) to the observa tional data. Drukier (1995) 
and Fryor et al. ( 1986a|b I have pointed out the limi- 
tations ot fitting globular clusters with single compo- 
nent models. We have tried to mitigate these issues to 
an extent. In particular the background potential in 
our case is supplied by an ensemble of lowered isother- 
mal distribution functions of different stellar masses in 
equilibrium such that the total number of stars in each 
mass bin is proportiona l to the observed mass func- 
tion ( |Richer et aL|[2008| >. As we shall argue, because 
our field is centered outside the half-light radius of the 
cluster, the effects of mass segregation are more mod- 



est than if we had observed the entire cluster or a field 
near the core, so a particular mass component (or value 
of a) is sufficient to characterize the data in our field. 
We use the best-fitting model to assist in deconvolving 
the projected distributions and providing an estimate 
for the mass of the cluster and the escape velocity for 
stars within our field. 

1.1. Observational Overview 

Two sets of observa tions with the Ad vanced Cam- 
era for Surveys (ACS, |Ford et al.l|1998[ ) on the Hub- 
ble Space Telescope (HST) of the globular cluster 
NGC 6397 spaced over five years provide sensitive 
probes of the dynamics of the globular cluster. The as- 
trometry is sufficiently sensitive to resolve the proper 
motion of individual stars in NGC 6397, probing the 
theoretical model of the cluster in detail as well as pro- 
viding constraints on its mass, distance and relaxation 
time. This stud y is closest in spirit to that of |McLaugh-| 



lin et al.| ( |2006[ ), so it is quite valuable to introduce our 
study by comparing and contrasting it with this work 
on stars in the core of 47 Tuc. Most importantly for 
this work, 47 Tuc has perhaps ten times more stars 
than NGC 6397 and the MAM06 sample focuses on 
the core of 47 Tuc, so MAM06 have about 13,000 stars 
in their proper motion sample while we have about 
3,000 stars. The larger number of stars in 47 Tuc 
and its larger physical size have the dynamical effect 
of giving the l arger cluster a long half- mass relaxation 
time of 4 Gyr (Gnedin k Ostriker 1997) compared with 
0.3 Gyr for NGC 6397 (§ |4.4|). The dynamical evolu- 
tion of 47 Tuc as a whole has not proceeded to the 
same extent as that of NGC 6397, so we expect to see 
different processes at work. Second, our sample lies be- 
yond the half-light radius of the cluster, so our study 
is more sensitive to the global properties of the cluster 
and allows us to make an estimate of its mass without 
large extrapolation. 

As we have argued, the work of MAM06 provides 
a natural benchmark and also a treasure map as 
well. Many of the figures, tables and arguments pre- 
sented here shall be reminiscent of MAM06; however, 
NGC 6397 and the pure ACS astrometry present many 
new challenges and new opportunities in the study of 
globular clusters. 

1.2. Outline 

Table [T] presents some of the fundamental properties 
of NGC 6397 as well as auxiliary information that af- 
fects our observations of this particular cluster. We 
have attempted to give the original references for the 
various values; however, as this is not intended to be a 
thorough review, this list may not be exhaustive. The 
table also contains results from this paper with the 
corresponding section and acts as an index of sorts to 
delve directly into the results of interest. When a re- 
sult relies on data presented here combined with other 
work (typically MM91), the second paper is also listed. 

The following section (§|2| outlines the obs ervations 



Next we describe the sample selection (§ 3.1 ), the field 
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TABLE 1 
Basic data on NGC 6397 



Property 



Value 



Reference 



= 17 h 40 m 42!09, S = 
338°1650,6 = -11' 
t = 5.73 



Cluster center (J2000) 
Galactic coordinates 
Apparent Magnitude 
Luminosity 
Integrated colors 
Metallicity 

Foreground reddening 
Foreground absorption 

Field contamination (V < 21) 
(V < 29) 
Structural Parameters: 

Total Mass* 

Core Radius 

Half-Light Radius 
Multi-Mass King Model: 

Central Escape Velocity 

Central Escape Velocity 

CT-paramctcr (19.5 < F814W < 24.5) 

(T-parameter (F814W < 16) 
Heliocentric distance: 

Sub dwarf fit 

White-dwarf fit 
Kinematic 
Timescales at 5': 
Crossing time from proper motions 
Evaporation time 
Relaxation time from white-dwarfs 
Dynamical relaxation time (in our field) 
Dynamical relaxation time (at R,^) 



-53°40'27'f6 
9595 



L v = 4.6 x 10 4 rf| 53 L 
B-V = 0.73,(7 - V = 0.85 
[Fe/H] = -2.03 ±0.05 
E(B - V) = 0.18 ±0.01 
A v = 0.56 
^4-F814W = 0.33 
Sf orc ft: 23 stars arcmin -2 
ft: 240 stars arcmin -2 

1.1 ± 0.1 X 10 5 d| 53 M Q 

3 arcseconds 

Rh = 2.9 arcminutes 

2.31 mas/yr 
2.8 <i^53 mas/yr 
1.01 mas/yr 
0.54 mas/yr 

2.53 ±0.05 kpc 
2.67 ± 0.25 kpc 
2.55 ±0.11 kpc 
2.0 ±0.2 kpc 



r/cr : 



Tc 

r e = (dlnN/dt) 
r r > 0.7 Gyr 
r r ft; 1 Gyr 
r rh ft: 0.3 Gyr 



0.6 Myr 

i ^ 



3 Gyr 



Goldsbury ct al 
Goldsbury ct al 
Trager et al.|l995 
'trager et al .l 



2010 



2011 



1995 



Harris|l996, 
Gratton etal.|2 0Q:; 



Reid & Gizis 1998 



Gratton et 



Cardelli et al.|l989 
Sirianni ct al. 2005 
Ratnatunga &: Bahcall 



aTj 



2003 



|1985 



Trager et al. 


1995 


'1'rager et al. 


1995 


This paper, 5 


4.1 


Gnedin et al. 


1999 


This paper, § 
This paper, § 


4.1 
4.1 



Gratton et al.|2003| 

Reid fc Gizis|l998l 

Hans eneral.]200_7] 

Meylan fc Mayor|l991 this paper, 



4.6 



This paper, 
This paper, 
This paper, 
This paper, 
This paper, 



All of the error intervals are ninety-percent confidence regions. 
2.53 kpc. 



*The value (I2.53 is the distance of NGC 6397 from Earth divided by 



geometry 
ticular § 3.3|out 



) and the analysis techniques. In par- 
ines the theoretical model and the sta- 



describes the esti- 
the results follows 



tistical estimators to probe it (an app endix applies this 
model to proper motions) , and § |3.4 
mation of er rors. A presentation ot 
m §|U § |4.1| presents a the oretical model for the current 
state of the cluster, § |4.2| ou tline s the evidence for mass 
segregation in the cluster, §|4.3| looks for anisotropy in 
the velocity distribution, § |4.4| describes the dispersion 
in proper motion as a fun ctio n of stellar type and posi- 
tion in the cluster 



and 



4.5 examines the proper mo- 
tion distribution in further detail for various subsam- 
ples. In each of these sections the results of the model 
are presented with the observations. We als o ob tain 
an esti mate of the distance of the cluster (§ 4.6 ), it s 
mass (§ 4.7 1 and the candidate stellar escapers [§ 4.8 1; 
all of these areas build directly upon the model of the 
cluster. Next, § [5] discusses the consequences of thes e 



results for our understanding of the d ynam ics (§ 5.1| 



the global properties of NGC 6397 (§ [5^ and fuTure 



directions (§ 5.3 1 



2. OBSERVATIONS 

The first epoch photometry used in this paper is the 
same as that in the various papers on this cluster we 



developed from our HST Cycle 13 project ([Ri c her et a~ 
2006| lAnderson et al.||2008| |Richer et al.||2008[ ). There 
are two ACS epochs that went into determining the 
stellar proper motions. Nine new orbits were added in 
2010 to re-image the cluster in F814W only, in order 
to obtain the best proper motions possible. Because 
all the astrometry was performed on F814W images 
and used the distortion corr ections for the F814W fil- 
ter ((Anderson & King 2006)), there should not be any 
color-dependent effects on the astrometry. These ob- 
servations are summarized in Tabled We constructed 
a sample of likely stars using the techniques outlined in 
Anderson et al.| ( |2008[ ) to remove artifacts and galaxies 
from the sample. 

2.1. Second- Epoch Astrometry 

We reduced the second epoch images using the pub- 
li cly available routine img; 2xym_WFC . 09x10, described 



. Anderson fc King] ([2006 ) . This is a one-pass photom- 
etry program that goes through each exposure pixel- 
by-pixel and identifies local maxima that are suffi- 
ciently bright and isolated. The routine then measures 
positions and fluxes for these putative stars by means 
of a spatially variable library PSF. We identified every 
single local maximum as a potential star and measured 
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it with the PSF. This gave us a list of 1.5M "sources" 
for each exposure. The vast majority of these are noise 
flucutations, but some are real stars. 

In order to determine proper motions, we must com- 
pare the positions of the stars in the first epoch with 
those in the second epoch. To do this, we must 
transform the positions of the stars in each expo- 
sure into the reference frame. We thus took a list 
of the main-sequence stars (i.e., those that lie along 
the main-sequence ridgeline, MSRL) and identified the 
same stars in each of the second epoch images. We 
took these pairs of positions to construct a general 
six-parameter linear transformation from the second- 
epoch exposure into the reference frame. We exam- 
ined the residuals of this transformation and noted that 
the distortion solution had changed somewhat over the 
time from the first epoch (in 2005 pre-SM4) to the sec- 
ond epoch (in 2010, post-SM4). This visible residual 
distortion could be characterized by clear linear trends 
(with an amplitude of ~0.03 pixel) within each chip, so 
we decided that instead of transforming the entire de- 
tecor with a global transformation, we would transform 
each chip indepndently with a global transformation. 
The residuals of this approach showed no significant 
remaining trends. We validate this with local transfor- 
mations below (§ 2.2 1. 

To compute the displacement for each star between 
the first and second epochs, we used the above trans- 
formations to map each "peak" detection in each im- 
age found above into the reference frame. We did this 
separately for the 252 F814W deep first-epoch expo- 
sures and for the 18 deep second-epoch exposures. To 
identify the most likely position of each of the 46,785 
sources in each epoch, we collected all the peaks found 
within five pixels of the target location and identified 
the place where we found the largest number of detec- 
tions within a radius of 0.75 pixel. We then performed 
an iterative sigma clipping (with a c-threshold of 3.5) 
to determine an average position and a error in the 
average, based on the RMS about the average and the 
number of peaks remaining. We did this for the first 
and second epochs. 

Although we started out with "good" positions for 
the first-epoch stars, the method we used to find the 
positions for these displacements was different enough 
from the original finding algorithm that we needed to 
be robust against the impact of artifacts and neighbors. 
For instance, a faint star within five pixels of a brigher 
star (or a PSF artifact) might have the position of the 
brighter star identified as its position, rather than its 
own position. To avoid this, we used the detection- 
concentration method to re-find the first-epoch posi- 
tions. If the new first-epoch position is not consistent 
with the original first-epoch position, then it is likely 
that the second-epoch position will also be bogus, so 
we flag these stars as not having valid proper motions. 
If, on the other hand, first-epoch data shows that the 
most significant concentraion of detections within five 
pixels is in fact the source itself, then we can be confi- 
dent that the second-epoch images will not mis-identify 
it as a neighbor or artifact, and we can trust the proper 



motion. 

To determine a displacement, we simply subtract the 
first-epoch position obtained above from the second- 
epoch position. The random error in the displace- 
ment is simply the sum in quadrature of the first- and 
second-epoch errors (dominated of course by the sec- 
ond epoch). The displacement errors for the bright 
stars are typically 0.002 pixel and those for the faintest 
detectable stars are 0.20 pixel. We convert this into 
a proper motion by dividing by the time baseline of 
4.936 years. We neglect the statistical error related 
to the global linear transformations, since they were 
based o n we ll over 1,000 stars in each chip. We show 
below (§ 2.2 1 that if we apply a local correction to the 
transformations, we get residuals that are at about the 
level of the quoted measurement errors. Much of this 
is due to the small number of stars used to do the local 
transformations and the internal dispersion. However, 
it is clear that transformation error is no larger than 
the quoted random measurement error. 

To measure the proper motions for the bright stars, 
we used the 40s short exposures from the two epochs. 
The first epoch has two 40s exposures and the second 
epoch has three. We reduced the pixel-based-CTE- 
corrected images as above for the deep exposures, then 
as before we found the transformations into the refer- 
ence frame using the common stars along the MSRL 
and global transformations for each chip. We then 
computed displacements from first to second epoch 
(and the corresponding errors) and report the proper 
motion in pixels per year. This adds 999 stars to the 
proper motion list. About 85 sub-giant-branch, red- 
giant-branch and horizontal-branch stars are saturated 
in the 40s exposures. Because of CTE concerns in the 
extremely low-background 5s and 10s exposures, we 
decided not to attempt to compute proper motions for 
these stars. 

We thus have an estimate of the proper motion for 
the faint stars from the set of long exposures and an 
estimate of the proper motions for most of the bright 
stars from the short exposures. This gives us motions 
for essentially all stars at or below the turnoff. Fig. [I] 
gives the proper motion relative to the mean cluster 
motion of all the stellar objects measured in the field 
with proper- motion errors less than 0.4 mas/yr. The 
arrow gives the direction of the cluster center from the 
ACS field center. It is coincidental that this arrow 
happens to point toward the field stars in the proper- 
motion diagram. 

2.2. Global vs. Local Transformations 

When we performed the single-chip global transfor- 
mations above, we did not notice any remaining trends 
of motion with position, indicating significant resid- 
ual distortion, but this is worth verifying. To do this, 
we followed the example of MAM06 and computed for 
each star a local correction for its proper motion based 
the average motion of the nearest 75 MSRL cluster 
stars. The typical cluster star has a random motion of 
0.095 pixel per year with respect to the systemic cluster 
motion, with an error of 0.0045 pixel per year (about 
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5%). We found that the median systematic proper mo- 
tion residual was 0.008 pixel per year. Since this was 
based on the average of 75 stars, we would expect to see 
0.007 pixel per year without any true systematic error 
being present at all. The fact that we observe slightly 
more than this could be indicative of non-random er- 
rors of 0.0045 pixel per year, similar to the random 
measurment error. At any rate, our choice of using 
global transformations to transform each chip into the 
reference frame will not have a significant impact on 
our motions. 

2.3. Charge- Transfer Efficiency Correction 

Before any positions were measured for the second 
epoch, the images were corrected for imperfect charge- 
transfer efficienc y (CTE) using the pixel-ba sed correc- 



TABLE 2 

ACS observations of NGC 6397 



tion described in Anderson & Bedin 
correction has been applied to the ca 



(2010J). The CTE 
jrated (-fit) 



ages, as is the practice in the ACS pipeline. The back- 
ground in the deep exposures is about 125 electrons 
per pixel. This large background shields faint sources 
from the majority of traps they would experience. The 
result of CTE losses would be to cause the stars to be 
shifted towards the gap, the faint stars more than the 
brighter stars. The single-chip-based transformations 
we used compensated somewhat for the average CTE, 
since it naturally allowed for an arbitrary rescaling of 
the y axis from first to second epoch. So all we are sen- 
sitive to is differential CTE shifts: a shift of the faint 
stars relative to the bright stars. 

To get a sense of the overall amplitude of the as- 
trometric CTE effect, we reduced the CTE-corrected 
and the CTE-uncorrected images in an identical way 
and compared the output positions as a function of 
the star's flux. We found that the brightest stars had 
a +0. 02-pixel relative shift at the gap, and the faintest 
stars had a —0.04-pixel shift. The CTE model is not 
perfect, but it ca n be counted on to remov e roughly 
75% of the signal (Anderson & Bedin 2010). As a re- 
sult, we are confident that the residual C l'E-related 
astrometric error will be less than 0.005 pixel in total 
displacement, which corresponds to a proper motion 
of 0.001 pixel per year. Such a small signature is not 
possible to see in the data, since the random motions 
of the stars are about 10 times this. 

2.4. Incompleteness Corrections 

In this paper we make no use of incompleteness cor- 
rections. We are not considering very faint stars in 



tions (Hansen et al. 2007 


Richer et al.|2008 


Anderson 


et al.|2008) the faintest stars in our current sample are 



epoch stars in our subsequent epochs; consequently, 
inco mpl eteness should not be a serious issue. In fact 



4.2 (Fig. fTJ) we present the completeness fraction 
as a function of magnitude to estimate how differential 
incompleteness could bias the radial distributions and 
find it to be much smaller than the statistical uncer- 
tainties. 



Data Set 


Program ID 


N ohs 


Filter 


Total Exposure 


Richer 2005 


10424 


126 


F606W 


93.442 ksec 


Richer 2005 


10424 


252 


F814W 


179.704 ksec 


Rich 2010 


11633 


18 


F814W 


24.705 ksec 




H a [mas/yr] 

Fig. 1. — Proper motions of the objects identified as stars in 
the ACS field. Since we used cluster members to define the 
reference frame, the zero point of the vector-point diagram cor- 
responds to the bulk motion of the cluster. Those stars within 
the inner red circle comprise the sample to determine the color- 
magnitude diagram of the cluster. The stars within the outer 
blue dashed circle comprise the sample used to probe the proper- 
motion distribution. The arrow indicates a proper motion in the 
direction of the cluster center from the center of the ACS field. 
The swarm of stars centered around fi a PS — 12 are field stars. 




[mas/yr] Total Motion 



F606W-F814W 



Fig. 2. — Proper motions, proper-motion errors and colors as 
a function of magnitude. The left panel give the observed total 
proper motions as a function of apparent magnitude in red, and 
the proper motion errors in blue. The right panel plots the 
apparent magnitude against the observed color. Fur thermore, 
the boxes depict the four subsamples outlined in §3.1| 
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3. METHODS 
3.1. Star Selection 

To construct a member-only color-magnitude dia- 
gram, we select only those stars with proper motions 
(Fig. IT and[2]) relative to the mean motion of the cluster 
stars of less than 2 mas/yr for 15.5 < F814W < 26 (red 
circle) and proper-motion errors (blue points in Fig. [2] 
and all subsequent figures) of less 0.4 mas/yr. All stars 
brighter than F814W = 15.5 are used to construct the 
color-magnitude diagram, but because proper motions 
were not determined for these stars, they do not end 
up in the proper-motion sample. For stars fainter than 
F814W = 26, the proper-motion cutoff is 5 mas/yr and 
there is no error cutoff, field stars outside of the clus- 
ter comprise the broader distribution centered around 
\n\ « 15 mas/yr on the left panel of Fig. [2j 

Selecting the stars moving with the cluster yields 
the color-magnitude diagram (right panel of Fig. pi 
from which the proper-motion sample is constructed. 
To prevent biasing the selected proper-motion sample, 
the cutoff in proper motion for this sample is 5 mas/yr 
and a proper- motion error cutoff of 0.4 mas/yr for 
all magnitudes. The sample is further winnowed by 
including only stars with colors that are within 0.05 
magnitudes of the main-sequence and the white-dwarf 
cooling track. Table [3] gives some statistical proper- 
ties of the objects that we identify as cluster and field 
stars. We focus on two large subsamples to probe the 
mass segregation of the cluster in radius. In partic- 
ular, sample Rl consists of main-sequence stars with 
17 < F814W < 19.5 that have a typical mass of about 
O.5M . The second sample R2 probes fainter MS stars 
(21 < F814W < 22.5) with a typical mass of O.2M . 
Here we want to maximize the size of the samples to be 
sensitive even to subtle differences in the radial distri- 
butions without allowing the two samples to overlap. 
The samples Rl and R2 represent a compromise. 

We also need to select some samples to investigate 
the dynamics of the cluster through the proper mo- 
tions. Our results concerning the proper motions of 
the population as a whole focus on those stars with 
the most precisely measured proper motions, the "Best 
PMs" sample, the main-sequence stars with 19.5 < 
F814W < 24.5 — all of these PMs were measured us- 
ing the long exposures, making a uniform sample. To 
probe the dynamical relaxation of the cluster we will 
define four smaller samples. Two of these are subsam- 
ples of the "Best PMs" sample, and two come from 
other regions of the color-magnitude diagram: 

1. White dwarfs with 22.5 < F814W < 25 ("Bright 
White Dwarfs"), 

2. Main-sequence stars with 22.5 < F814W < 25 
("Faint MS"), 

3. Main-sequence stars with 19.5 < F814W < 20 
("Middle MS"), 

4. Main-sequence stars with 15.5 < F814W < 17.38 
("Bright MS"). 



Spherical Radius 
Cylindrical Radius 




Radius from Center farcminutes 



Fig. 3. — The projected (blue dashed) and spherical radii (red 
solid) probed by the ACS field. 

Sample 1 contains 46 white dwarfs with the best mea- 
sured proper motions. Sample 2 contains 254 main- 
sequence stars with well-measured proper motions over 
the same range of apparent magnitude as the white 
dwarfs. The size of sample 2 sets the goal size for the 
remaining subsamples. Sample 3 consists of 255 main- 
sequence stars, the brightest stars with proper motions 
measured from the long exposures. Furthermore, these 
main-sequence stars are expected to have a mass sim- 
ilar to the white dwarfs. Finally, sample 4 consists 
of the brightest 255 main-sequence stars with proper 
motions (here measured in the short exposures). This 
final sample provides our best estimate of the dynam- 
ical properties of the progenitors of the white dwarfs. 
The three main-sequence samples have the nearly the 
same size, allowing a fair comparision each sample with 
the white-dwarf sample. All four subsamples are de- 
picted in the color-magnitude diagram. Table [3] gives 
the astrometric errors of the entire sample as well as 
the smaller subsamples. 

3.2. The ACS Field 



Fig 



To efficiently average over the ACS field (depicted in 



20 



with the potential stellar escapers), we define 
a window function that quantifies the regions of the 
cluster probed by the field. In the plane of the sky, the 
window function is simply defined to be proportional 
to the ratio of the area of the ACS field at a given pro- 
jected radius from the center to the total area at that 
radius. The geometry of the field provides the initial 
guess of the window function. The estimate of the win- 
dow function is optimized using the assumption that 
those stars in the field that are not cluster members 
are uniformly distributed. 

For convenience the window function is normalized 
to unity at its peak. The three-dimensional window 
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TABLE 3 

Proper-motion dispersions, root-mean-squared errors, mean errors for the various subsamples 







Median 






Mean Error 


RMS Error 


/ 


Sample 


Number 


F81 4W 


Mass 


[mas yr — 1 ] 


[mas yr — 


[mas yr ] 


[mas yr" 1 ] 


[%] 


All stars 


6345 


21.6 




5.42 


1.31 


0.07 


0.10 


0.30 


All MS stars 


2880 


20.3 


0.34 


0.36 


0.37 


0.03 


0.04 


0.58 


All WD stars 


186 


26.1 


0.53 


0.50 


0.47 


0.23 


0.25 


15.67 


All field stars 


3279 


23.3 




4.40 


2.79 


0.09 


0.12 


0.10 


Radal Samples: 


















Rl (MS: 17 < F814W < 19.5) 


788 


18.3 


0.58 


0.37 


0.37 


0.05 


0.05 


1.09 


R2 (MS: 21 < F814W < 22.5) 


736 


21.6 


0.18 


0.38 


0.37 


0.03 


0.03 


0.28 


Proper-Motion Samples: 


















Best PMs (MS: 19.5 - 24.5) 


1899 


21.0 


0.24 


0.37 


0.37 


0.03 


0.03 


0.30 


Bright WD stars (1, 22.5 - 25) 


46 


24.3 


0.53 


0.34 


0.32 


0.08 


0.09 


4.15 


Faint MS stars (2, 22.5 - 25) 


254 


23.1 


0.11 


0.37 


0.38 


0.05 


0.05 


0.83 


Middle MS stars (3, 19.5 - 20) 


255 


19.8 


0.42 


0.34 


0.35 


0.02 


0.02 


0.23 


Bright MS stars (4, 15.5 - 17.38) 


255 


16.7 


0.74 


0.34 


0.34 


0.04 


0.05 


1.04 



The quantity / = 1 — Sv^true/on bs gives the relative decrease in the value of the velocity dispersion after correcting for the uncertainty 
in the observed proper motions. The "All Stars" sample has no proper-motion cutoff; whereas the non-fieldx other samples have a 
proper-motion cutoff of 5 mas/yr. All samples have an proper-mot ion-error cutoff of 0.4 mas/yr. The masses of main-sequence stars 
are given in solar masses using the models oflHurle y et al . (2008) and assuming a distance of 2.53 kpc. The masses for the white 
dwarfs come from the spectroscopic measurements ol |rvalirai eT al. (2008). 



function is defined similarly: the ratio of the total vol- 
ume within the field as extended through the cluster 
at a given three-dimensional radius from the center to 
the total volume at that particular radius. Again this 
function is normalized to unity at its peak. These win- 
dow functions, as depicted in Fig. |3j provide a useful 
technique to generate either a Monte Carlo sample or 
integrate over the three-dimensional structure of the 
cluster. The window function is simply the probability 
that a particular star in the entire Monte Carlo sample 
ends up in the sample that represents the field. 



3.3. Theoretical Model 

To obtain estimates of the total mass of the cluster 
and the escape velocities of individual stars, we will 
rely on a theoretical model of the cluster. We will use 
the model as one method of deprojecting our velocity 
and surface density measurements to constrain the dy- 
namical properties of the cluster. This section outlines 
the model and the statistical measures that we will 
used to characterize both the model and the data. An 
appendix derives several useful results from the model 
that further our interpretation of the data. 

To model the stars in the globular cluster we will as- 
sume that the phase-spac e density is given by a low ered 
Maxwellian distribution ( |Michie||1963| |King||1966[ ), 



f(E) = N [exp(~E/<T 2 ) - 1] 



(1) 



where N is the normalization, E is the energy of a 
star per unit mass (less than zero because the stars are 
bound) and a is a parameter with dimensions of veloc- 
ity that sets the velocity scale for a particular group 
of stars in the cluster. Let \I/(r) be the negative of 
the potential energy, which vanishes at the edge of the 
cluster (the tidal radius, r t ) so the number density of 



stars at a radius r is 
Pi 



n(r) 



(27TCT 2 ) 



3/2 



4ttv 



dv 



(2) 

where p\ is the number density normalization, v is the 
speed of the star and v 2 = 2^>(r ) is the speed needed to 
escape the cluster by reaching the tidal radius. There- 
fore, the number density is given by 



n(r) = pi 



V2c 



2v R 



1 



\_v_i 
3ct 5 



( 3 ) 

Similarly the local three-dimensional velocity disper- 
sion is given by 

V\r) = ^ (4) 
Jo 



where 



M-v 2 )/(2a 2 ) _ 1 



,,"+2 



dv. 



(5) 



To build a model star cluster we take a range of stel- 
lar masses each with in its own value of a in Eq. [3j 
the stars of different masses are assumed to be in ther- 
mal equilibrium, i.e. n%i<jf is a constant. This pro- 
vides the total mass density of stars as a function of 
the local potential. We select a value for the central 
potential and integrate the Poisson equation for the 
gravitational potential using the total mass density of 
the ensemble of stars as a source. For each value of 
the central potential, we vary the normalizations of 
the mass densities (pi) for each ma ss to match the ob - 
served mass function of NGC 6397 ( |Richer et alj2 008). 
This generates the value of the potential as a function 
of radius, yielding the number density, column density 
and the one-dimensional velocity dispersion as a func- 
tion of radius for any subpopulation of stars with a 
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particular value of a. We then choose the values of 
the central potential and a that yield the best match 



to the photometry of NGC 6397( |Trager et al.l|1995 |. 
§4.1 1 compares the observational results from this work 
and previous work with the theoretical model of the 
cluster outlined here, and the appendix derives several 
properties of this model that are useful to interpret the 
observational results. 

We would like to compare these model distributions 
against the data and characterize these distributions 
using standard deviation of the velocity distribution 
along a particular direction. The standard deviation 
has several advantages including the ease of calculation 
from both the models and the data and a straightfor- 
ward physical interpretation, in the Jeans equation, for 
example. However, the standard deviation is extremely 
sensitive to outliers. In fact a single interloping star, 
say from the Galaxy, in the sample could ruin our esti- 
mate of the standard deviation of the stars within the 
cluster. There are several robust estimators of the scale 
or width of a distribution, such as the median absolute 
deviation or for two-dimensional distributions like the 
proper motion, the median magnitude of the proper 
motion. These two estimators are maximally insensi- 
tive to outliers; up to half of the points can be shifted 
to arbitrarily large values without affecting the result. 
For a normal distribution these estimators yield values 
proportional to the standard deviation; however, as 
the distributions deviate from normality such as those 
present in the lowered isothermal sphere, these con- 
stants of proportionality can change. 

Instead of these estimators we use the first quartile 
of the differences of the proper motions along a partic- 
ular direction to estimate the standard devi ation of the 
distribution: (t he estimator Q n defined by Rousseeuw 
fc Croux|p93 l 



a = d n x first quartile of — fij\ : i < j) . (6) 

where fii and fij are the measured proper motions of 
stars in our sample or the model along a particular di- 
rection and d n is a factor that depends on the size of the 
sample that ensures the a is an unbiased estimator of 
the standard deviation for normally distributed data. 
As the sample size diverges (n —¥ oo), d n approaches 
l/(2erf _1 (l/4)) « 2.219. Like the medians discussed 
earlier, the estimator a has a breakdown point of 50%, 
meaning that one could shift up to one half of the val- 
ues Hi to infinity without affecting the estimator. 

The estimator a has several further advantages. 
First, it can be determined for each direction, unlike 
the median magnitude of the proper motion. Sec- 
ond, the value of a approximates well the standard 
deviation for distributions that differ significantly from 
Gaussian such as a uniform distribution. Finally, the 
estimator a is statistically efficient, meaning for a sam- 
ple of a given size the typical error is smaller than that 
which results for the median absolute deviation or me- 
dian proper-motion magnitude. 

We calculate <x from the proper-motion data, the ra- 
dial velocity sample and from the models. Additionally 
we treat it as a proxy for the standard deviation for 



much of our analysis. In particular, when we correct 
our observed proper motions for the estimated errors 
in the proper motions, we will use 



~ 2 _ » 2 

^/i,truc ^^i.obs 



1 " 

n ^ 

i=l 



(7) 



where is the proper-motion error of each star. Ta- 
ble [3] shows that these cor recti ons are generally small. 
Finally as we shall see in § |4.7| the Jeans equation uses 
the velocity variance to provide an estimate of the mass 
of the cluster; therefore, we shall use both observed 
proper-motion variances and a to estimate the mass of 
NGC 6397. We shall find that the values of a and the 
standard deviation are approximately equal for most 
of the subsamples that we investigate, so this choice 
makes little difference to the final mass estimates. 

3.4. Error Estimation 

The statistical bootstrap, the cousin of the 
Qucnouil le-Tukey jac kknife, was introduced by Efron 
in 1977 (Efron] 1979). In this section we introduce or 
rather reintroduce the bootstrap and outline how we 
use it to estimate the errors in our derived proper- 
motion dispersions, colum n density di stributions and 
cluster mass. (Refer to |Lupton||1993| for further de- 
tails) . Observations of NGC 6397 draw a sample of n 
stars with positions and velocities {(^,Wi)}o from the 
phase-space distribution function, f(x,v 7 to). We use 
various statistics of the observed stars to make conclu- 
sions about the distribution function and its evolution: 
for example, the mean position along the a-axis of the 
sample: 

1 x - 



and its associated error 
1 1 



-1^ 



(8) 



(9) 



By performing these summations, we assumed that the 
statistics of the set {(Hi, Vi)} are a good approximation 
to those of the actual phase-space distribution func- 
tion, f(x,v,t ). This is also the assumption behind 
statistical bootstrapping. 

However, we could have taken a more complicated 
route and created another distribution function (in fact 
the distribution that maximizes the likelihood of the 
observations) , 



Xi,V- Vi), 



(10) 



and then calculated mean position by taking the inte- 
gral over the region probed by the ACS observations, 

x a f*dxdv. (11) 

'ACS Field 

Now we have two methods to estimate the error in 
our determination of the center of mass: calculating 
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an integral analogous to the second summation (Equa- 
tion^]), or drawing a variety of samples from our new 
distribution function /* and investigating the distri- 
bution of x* a . Our original sample, {(xi,Vi)}, is now 
only one of many possible realizations of /*, and the 
distribution of x* over the various resam plings approx- 
imates the distribution of x* a over /* ( |Lupton| [l993) 
which we in turn assume to approximate / (a; , v, t f ). 
Thus, through resampling we can estimate the error 
in our determination of 5* given by Equation [9] This 
Monte-Carlo resampling is less expedient than using 
Equation [91 but often we do not have the luxury of a 
statistic denned as simply as x a is in Equation|H] Even 
more rarely do we have the luxury of a straightforward 
analytic error estimator like Equation [9] 

So how is the bootstrap implemented with ACS ob- 
servations? Ordinarily when one calculates a physical 
value from a sample, one uses all the stars, counting 
each one once. To use the bootstrap, one simply needs 
to select from an n-star sample, n stars with replace- 
ment. To create a bootstrapped sample, we uniquely 
label each of the stars in an n-star sample from 1 to 
n. We then select n integers from 1 through n and use 
the stars with the appropriate labels to create a new 
ensemble. Some of the stars in the original sample are 
included in the new ensemble several times and others 
not at all. We often use the value of a for a sample 
of proper motions to characterize them. We resamplc 
each group of stars 1,000 times, determine the value 
of a of each resampling, sort the results and quote the 
fifth and ninety-fifth pe rcen tiles as a ninety percent 
confidence region. In § |4.7| we estimate the mass of 
NGC 6397. In this case we calculate a mass estimate 
for each resampling, sort these results and calculate the 
ninety-percent confidence region for the final result. 

4. RESULTS 
This sec tion begins with a theoretical model for the 



cluster (§ 4.1 1 to provide a context for various char- 



acteristics ofthe data inclu ding the radial distribution 



(§ |4.2[), proper motions (§§ 4.3 4.5 1 and some der ived 



properties o f the se data including its distance (§ 4.6) 
and mass (§ 4.7| 



4.1. A Model for NGC 6397 



The mass function of 
served concentration O] 



Richer et al.|(|2008) and the ob- 
the cluster from Trage r et al. 



(1995) provide a useful starting point to construct a 
preliminary multi-mass King model for the cluster from 
which we can draw some global inferences. Fig. [4] com- 
pares the model surface density for the best fitting val- 
ues of the central escape velocity and dispers ion pa- 
ramet e r (a) with the observed photometry from Trager 
et al. (1995). W e ha ve also included the star counts 
within our field; § |4.4| and Fig. [13] presents these counts 
in further detail. To convert the star counts to an es- 
timate of the surface brightness we convert the ACS 
magnitudes used here to the Johnson system usin g 
the synthetic transformations of |Sirianni et al. (2005). 
This yields the total flux from ail ot the stars in the 
CMD sample defined in Fig. [2] of F tota i = 9.75 within 



the ACS field of 10.9 square arcminutes. This gives 
a mean surface brightness of 21.23 V magnitudes per 
square arcsecond or equivalently a mean V— magnitude 
of 17.94 for each main-sequence star in the proper- 
motion sample used to calculate the results depicted in 
Fig. 13 The star counts as a function of projected ra- 



dius are scaled by this mean magnitude per star, yield- 
ing the blue squares in Fig. [4] 

The radial, magnitude, velocity normalizations, cen- 
tral velocity dispersion and a— parameter of the models 
are arbitrary, and a comparison with the data is u sed to 
set the scale of t he model. We use the models of IHur-l 
ley et al. (20081 to set the ratio of the a— parameter 
ot the turn-oft stars, which are assumed to provide the 
bulk of the light, and that of the proper-motion sample. 
We take the mass of the turn-off stars to be 0.84 M Q 
and the proper-motion sample to be 0.24 M . Fur- 
thermore, we assume that the values of a following 
the equiparition relation with same value rriiaf to the 
various populations. The values of a and /z 6i o are de- 
termined with an unweighted least-squares fit between 
the models and the dat a. The bulk of the po ints in this 



fit are provided by the Trager et al. (11995) photome 



try; however, a measurement ot the velocity dispersion 
is required to set the velocit y sca le of the models. Here, 
we also use the data from § |4.4| and Fig. [13] Fig. [5] de- 
picts the observed proper motions, thosefound in the 
model and the observed radial velocities of MM91 for 
different assumed distances. 

Both the gross properties of the cluster such as its 
mass and the profile near the half-light radius depend 
little on the properties of the central regions. Op- 
erationally the radial scale and the ratio of the cen- 
tral escape velocity to the er-parameter are determined 
mainly by the surface photometry (Fig. 4| and the ve- 
locity normalization is determined by the value of a 
as a function of radius (Fig. 5|. The two crucial pa- 
rameters for the model are the central escape proper 
motion of /i ei o — 2.31 mas/yr and the er-parameter 
of 0.54 mas/yr for the turn-off stars that dominate 
the light and 1.01 mas/yr for the stars in the proper- 
motion sample. For the brightest stars, the value of 
the cr-parameter is approximately the central velocity 
dispersion along one dimension because [i e ^ a. The 
ratio of these two quantities is determined mainly by 
the surface brightness profile while the normalization 
is determined exclusively by the proper-motion data 
depicted in Fig. [5] The observed velocity dispersion, 
either in the radial direction or in a particular direc- 
tion of proper motion, is integrated over the line-of- 
sight so we can project the model onto the plane of the 
sky. Figure [5] depicts both of these results. Beyond 
one half-light radius the observed velocity dispersion 
depends only very weakly on the cr-parameter, so al- 
though mass segregation is important in generating the 
radial distribution, the velocity dispersions are nearly 
independent of stellar mass. Furthermore, the escape 
velocity at a given radial distance from the center of 
the cluster is typically three times the velocity disper- 
sion within five-percent agree ment with the analytic 
treatment in the appendix (Eq. A7). 
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Projected Model 
fi e ^ — 2.31 mas/yr 

E (cr — 0.54 mas/yr) 

£ (a — 1.01 mas/yr) 

Data 

Trager et al. 
■ This Paper 



Projected Radius [arcmin] 



FlG . 4. — Surface photometry as measured by |Trager et al.| 
p995|l (cyan crosses) and surface star counts from this work 
(blue squares) compared the the model surface density for the 
best fitting stars in the model (red solid curve) and the most 
massive stars in the model (blue dashed curve). 



Projected Model 
/i ej o — 2.31 mas/yr 

£r„ (a — 0.54 mas/yr) 

ct^ (a — 1.01 mas/yr) 

M./3 

Data 

This Paper 
« MM91 at 2.53 kpc 
MM91 at 2.0 kpc 




Projected Radius [arcminutes] 

FlG. 5. — Projected velocity dispersion as a function of pro- 
jected radius: model (best-fitting stars in red and most massive 
stars in green), radial-veloci ty and proper-motion meas urements. 
The velocity measured of |Meylan fc Mayor| ( |1991} (MM91) 
have bee n converted to prop er motions assuming a distance of 
2.55 kpc | |Richer et al.|2008[ l and 2 kpc to bring agreement be- 
tween the proper-motion and radial-velocity measurements. The 
error bars from this paper are ninety-percent confidence inter- 
vals obtained by bootstrapping. The errorbars on the MM91 
data are 1.64-tr. 



Fig. [4] shows that the change in the slope of the sur- 
face brightness profile as a function of the er-paramcter 
(or mass) is largest around one arcminute from the cen- 
ter of the cluster. The difference in slopes is more mod- 
est at the projected radius of the ACS field; therefore, 
the mass segregation in the model is weaker here. The 
analytic treatment outlined in the appendix although 
illustrative, is not strictly applicable here, because the 
escape proper motion is approximately 1 mas/yr and 
the cr-parameter is 1.01 mas/yr, so the value of a ap- 
proximates that of \i e and does not greatly exceed it. 
The field is not in the regime where the distribution 
function is strictly linear, however, as we can see from 
Fig-0 

The model velocity and distance normalizations also 
yield a model mass for the cluster of 



M = 1.1 x 10 5 <^ 53 M 



(12) 



where 1^2.53 is the distance to NGC 6397 from Earth 
divided by 2.53 kpc. We can combine the mass esti- 
mate in Eq. 12 with the mean stellar mass of about 
to* = 0.35 Mo( |Dotter et al.||2007| jRicher et al.||2008 ) 
to obtain an esti mate of the relaxati on time at the 
half-mass radius ( Spitzer fc Hart|[l971 ) 



r rh = 0.138- 



M 1 / 2 R 



.3/2 



K GV2 i n (0.4M/m*) 



(13) 



to yield T r h — 0.45g?2 53 Gyr where we have sup- 
pressed the additional logarithmic dependence on dis- 
tance through the Coulomb logarithm. 

4.2. Mass Segregation 

Fig. [6] depicts the cumulative radial distribution of 
several groups of stars in the sample, including two of 
the subsamples described in § |3.1| and Table [3j If we 
compare the main sequence stars with 17 < F814W < 
19.5 and 21 < F814W < 22.5, we obtain a p-value 
of 7 x 10~ 5 by the Wilcoxon rank-sum test and twice 
this value for the Kolmogorov-Smirnov (KS) test, indi- 
cating that it is highly unlikely that these two groups 
of stars follow the same radial distribution. The two 
main sequence samples show significant segregation, 
even greater than the model distributions. The model 
distributions are typically more centrally concentrated 
than the stars in the ACS field because t hey a re based 
on the surface photometry as outlined in § |4T| Because 
the field lies ou tsid e the half-light radius, following the 
arguments of § |3.3| we do not expect mass segregation 
to be strong in this region and the model exhibits only 
weak mass segregation. Either the cluster has stronger 
mass segregation than the model, or perhaps the sam- 
ples are incomplete especially for fainter stars near the 
center. Because the incompleteness fraction is low (less 
than 7% for the stars in these two samples), we can dis- 
count the second possibility. Although we will present 
the observed column density of stars within the field 
and use this to estimate the radial density gradient, 
we use this gradient as input for only two of several in- 
dependent mass determinations in § |4.7| the results of 
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Models 

Bright White Dwarfs 
Faint Main Seqnenee 
17 < F814W < 19.5 
21 < F814W < 22.5 




— t- QJ 



Projected Radius from Center [arcniinutcs] 



F814W 



Fig. 6. — Cumulative projected radial distribution for various 
subsamples. The three model curves give the best-fitting model 
in the middle (a mass of 0.35 Mq or F814W ~ 21) and models for 
two other masses assuming that rrt;<7? is constant. The models 
of potter et al.| < |2007^ give the correspondence between masses 
and magnitudes. The left-hand model traces the distribution for 
a stellar mass of 0.58 M (F814W ~ 18), and the right hand for 
a mass of 0.18 Mq (F814W ~ 22). 



Fig. 7. — Mean projected radius of stars as a function of 
F814W. The open circles give the mean radius for the input 
artificial stars, and the closed circles give the mean radius of 
those recovered. The errorbars here and in subsequent figures 
represent ninety-percent confidence limits and are obtained by 
bootstrapping. The solid curve gives the completeness fraction 
as a function of apparent magnitude as calculated from the ar- 
tificial star tests. This curve is read using the right-hand axis. 



which agree with the mass determination from surface 
brightness measurements. 

Fig. [6] serves a second purpose. Although the faint 
MS and bright WD samples span the same magnitude 
range, the stars in the WD sample are typically fai nter 
than in the MS sample as shown in Table [3] In § |4.5| 
we shall argue that these two samples have significantly 
different proper-motion distributions; in principle, this 
could result from the two samples having different ra- 
dial distributions. The typically fainter white dwarfs 
might be more difficult to find toward the center of 
the cluster, biasing this sample outward and biasing 
the proper-motion distribution toward smaller values. 
The cumulative distributions depicted in Fig. [6] demon- 
strate that this is unlikely. A comparison of the radial 
distribution of the two samples yields a p— value of 89% 
by the KS test and 92% by the Wilcoxon rank-sum test. 

We also compare the radial distribution of the 46 
bright white dwarfs with the two larger main sequence 
samples and find probabilities of 98% and 48% from the 
KS tests, and 92% and 32% Wilcoxon rank-sum prob- 
abilities. Furthermore, when we compare the radial 
distribution of bright white dwarfs (22.5 < F814W < 
25) with a sample of 97 fainter white dwarfs (25 < 
F814W < 26.5 - not shown in Fig. f§, we find 54% 
KS and 51% Wilcoxon probabilities. This fails to sup- 



port the conclusions of Davis et al. ( 2008b ) which used 
a smaller sample to conclude that the distribution of 
young white dwarfs was radially extended relative to 
that of older ones and argue that a possible explanation 



for this phenomenon was a white-dwarf natal kick. 

A possible explanation for the differing distributions 
as a function of apparent magnitude is differential in- 
completeness. It is more difficult to find faint stars 
in more crowded areas of the sky near the centre of 
the cluster. We have addressed this question in sev- 
eral ways in Fig. [7j First we have looked at the stars 
in the image and sorted them by apparent magnitude 
into bins of 50 white-dwarf stars or 200 main-sequence 
or field stars. A priori we do not expect the field stars 
(at any magnitude) to be centrally concentrated with 
respect to the cluster, so a change in the mean radius 
of the field stars would have to result from differential 
incompleteness. Because the number of field stars is 
limited to about 3,000, to get a better statistical han- 
dle on this effect, we inserted nearly 200,000 artificial 
stars throughout the image of various magnitudes and 
tried to recover them. The recovery fraction depicted 
by the solid curve in Fig. [7] gives an estimate of the 
completeness fraction. The mean radius of the input 
stars as a function of magnitude is given by the open 
circles. The mean radius of the recovered stars is given 
by the closed circles. We have sorted the artificial stars 
into bins of 5,000 objects. We see that fainter than 
F814W « 23 the recovered stars are on average fur- 
ther from the cluster centre than the input stars, the 
hallmark of differential incompleteness by about 0.05 
arcminutes; however, the effect is much more subtle 
than the change in the mean radius along the main 
sequence and white-dwarf tracks. It is unlikely to be 
strong enough to account to the differences in the ra- 
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Angle between Proper Motion and Direction to Centre [Degrees] 



Fig. 8. — Histogram of proper-motion angle relative to the 
center of NGC 6397. 

dial distirbutions depicted in Fig. [6] where the median 
radius of the models differs from the data by more than 
a tenth of an arcminute.. 

4.3. Proper-Motion Isotropy 

We measure the angle of the proper motion rela- 
tive to the cluster center of the main-sequence stars 
with 19.5 < F814W < 24.5, proper-motion errors less 
than 0.4 mas/yr and total proper motions less than 5 
mas/yr. The distribution of this angle is depicted in 
Fig. IS] There is no evidence for anisotropy from this 
distribution — a KS test indicates that the cumulative 
distribution of 35% of samples drawn from a uniform 
distribution will differ from uniform at least as much 
as the distri bution in Fig. [SJdoes. The p— value for the 
Kuiper test ([I960, a generalization of the KS test for 
a distribution on a circle) is smaller at 13%. If we look 
at the direction of the proper motions relative to right 
ascension and declination, the p— value for the KS test 
increases to 39% and for the Kuiper test to 32%. The 
directions alone do not reveal any anisotropy. 

We can also compare the velocity dispersion along 
the radial direction (<tr) to the tangential direction 
(<7t) for the main-sequence stars with the most pre- 
cisely measured proper motions (19.5 < F814W < 
24.5) and find that a (radial) /a (tangential) = 0.96 ± 
0.04, i.e. just consistent with isotropy at 90% confi- 
dence. To examine this in greater detail we sort the 
main-sequence stars and white-dwarf stars by appar- 
ent magnitude into subsamples of 200 stars for the 
main sequence and 50 stars for the white dwarfs (the 
faintest subsample may contain fewer). The proper 
motions are measured in the radial and tangential di- 
rections (relative to the center of the cluster) and the 
ratio of the dispersion in the radial direction to that 
in the tangential direction is determined and depicted 
in Fig. |9j We can divide the sample by the projected 
radius (Fig 



10) or total proper motion (Fig. 11). As 



for the distribution of all the proper motions for stars 
with 19.5 < F814W < 24.5, we see no strong evidence 
for anisotropy in any of these subsamples. We also 
relaxed the assumption that the proper-motion ellipse 
is aligned with the radial or tangential direction and 
even in this more general case found no evidence for 



I " 1 1 1 " I " 1 1 1 " I " 1 1 " 1 1 " 1 1 



1 1 " I 



1 1 1 1 1 1 



Main Sequence 
White Dwarfs 



F814W 

Fig. 9. — The radial anisotropy of the proper-motion ellipse as 
a function of apparent magnitude for main-sequence and white- 
dwarf stars. Only stars with total proper-motion error less than 
0.4 mas/yr are included. 

anisotropy in either the entire sample or the subsam- 
ples. The aforementioned results describe the proper- 
motion ellipse relative to the center of the cluster. We 
can also examine the proper-motion ellipse relative to 
the field boundaries. For the total proper-motion sam- 
ple and the various subsamples, we do not find any 
evidence for anisotropy of the proper motions relative 
to the boundary of the field. 

4.4. Proper-Motion Dispersion 

To examine the proper-motion dispersion, we sort 
the main-sequence and white-dwarf samples by appar- 
ent magnitude into subsamples of 200 and 50 stars, as 
before. We calculate the value of a for each subsample 
in the radial (<t_r) and tangential directions (ctt) an d 
define the one-dimensional value of 



(14) 



Because the proper motions are nearly isotropic, the 
value in either direction is well approximated by the 
one dimensional value of a as defined above. Further- 
more, for a spherically symmetric stellar system even if 
the velocity distribution is not isotropic , by extending 
the results of Leonard k. Merritt] ( |1989[ ) we find that 

a% + 4 = 2a z 2 (15) 

where a z is the velocity dispersion along the line of 
sight; therefore, Eq. [14] gives the natural value to 
compare with radial-velocity measurements, such as 
MM91. The root-mean-squared error in the proper 
motion is scaled and subtracted from each estimate of 
the dispersi on ( and its confidence interval) in quadra- 
ture (see § 3.3 for more details). Fig. 12 shows that 
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Fig. 10. — The radial anisotropy of the proper-motion ellipse as 
a function of projected radius from the center for main-sequence 
and white-dwarf stars. Only stars with total proper-motion error 
less than 0.4 mas/yr are included. 
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Fig. 12. — The value of a as a function of apparent magnitude. 
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Fig. 11. — The radial anisotropy of the proper-motion ellipse as 
a function of total proper motion for main-sequence and white- 
dwarf stars. Only stars with total proper-motion error less than 
0.4 mas/yr are included. 

the proper- motion dispersion is nearly constant, pos- 
sibly increasing slowly with fainter magnitudes at the 
bright end and nearly constant for F814W > 19 along 
the main sequence. Such a weak trend is expected for 
a sample beyond the half-light radius of the cluster be- 
cause the velocity distribution is cut off by escaping 
stars as shown by the model curve. The trend in the 



data is a bit stronger than the model. The observed 
proper motions of the brightest fifty white dwarfs are 
smaller by nearly two standard errors than those of the 
faintest eighty main-sequence stars — <Jwd — &ms — 
—67 /ias/yr with a standard error of 36 /ias/yr. The 
distribution of the differences in the dispersion over 
the bootstrapped resamplings is approximately nor- 
mal and the probablility that &wd > &ms is about 
three percent. The main-sequence stars with masses 
similar to the white dwarfs have F814W « 18.8. If 
we take the two hundred main-sequence stars with a 
median magnitude of 18.8 and calculate <jwd — fc'MS! 
we obtain —90 /ias/yr with a one-percent chance that 
&wd > &MS- If we take the brightest two hundred 
main-sequence stars, we obtain a closer agreement with 
o'wd — & ms = — 26 /ias/yr and a 17% chance that 
&wd > &MS- If we take the brightest fifty main- 
sequence stars, the difference decreases to —14 /ias/yr, 
and the probability increases to 38%. Looking along 
the main sequence to brighter magnitudes, it is appar- 
ent that the best agreement is with the brightest stars 
in the proper-motion sample; these stars most resemble 
the progenitors of the white dwarfs. 

With a sample of about 2000 main-sequence stars 
with 19.5 < F814W < 24.5, we sort the stars in pro- 
jected radial distance from the center of the cluster 
and determine the dispersion as a function of radius as 
shown in Fig. |13| In this case each subsample consists 
of 200 stars. Again the confidence intervals are deter- 
mined by resampling, and the root-mean-square error 
i n pro per motion is subtracted from a as described in 
§ |3.3| Both the density and proper-motion dispersion 
decrease with increasing radius. The curves give the 
best-fit power-law functions of radius for the value of 
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Fig. 13. — The value of cr and stellar density as a function 
of projected radius. The red and purple curves trace the best- 
fitting power-law relations in radius to the data of the same 
color. The upper green curve traces the value of the proper- 
motion dispersion as a function of radius according to the model. 
The lower curve depicts the value of the projected stellar density 
according to the model. 

<7 and column density. We have 

,-0.21±0.09 



a « (0.37 ±0.06) 



R 



(0.048 ±0.030) - 
* 5' 



mas yr, 



D \-1.40±0.38 , 

R \ stars 



arcsecond 2 



(16) 



We can also use the proper motions and positions 
themselves without binning in radius to determine the 
best fitting power-law relations for the column density 
and velocity dispersion, yielding 



^ 0.36 



R 



mas yr, 



-1.32 

— J stars arcsecond. -2 (17) 

in good agreement with the fit ted parameters for a 
and the column density in Eq. ( 16 ) and the velocity 
dispersion fitted from Fig. |j (e^ 0.36(i?/5')~ ' 20 ) 

For clarity we have not quoted errors on these pa- 
rameters because they are intermediate results in the 
det ermi nation of the mass of NGC 6397, which appears 
in § |4.7| with error bars. Although the power-law fit to 
the value of a characterizes the data well, the fit for 
the projected stellar density is poor. The observed star 
counts vary up and down with radius where any rea- 
sonable model would predict a monotonic decrease in 
radius. This behaviour appears to be larger than ex- 
pected given the errors on the densities, so it probably 
reflects either inhomogeneities in the cluster or under- 



estimated error s in our analysis. However, the models 
developed in § |4.1| rely only weakly on the projected 
stellar densities across the ACS field because they also 
account for the surface-bri ghtne ss profile of the entire 



cluster from Trager et al. ( 1995 ) . 

The stellar density and velocity dispersion can pro- 
vide a n estimate of the rela xation time for a group of 
stars, ( [Spitzer fc Hart]|1971| ) 



0.065t;| MS 



3G 2 m>ln(0.4M/m») 



(18) 



where urms is root-mean-square velocity of the stars, 
n is the number density, M is the total mass of the 
cluster and to* is the mass of a typical s tar (taken to 
be 0.35 M ). Using the results of Eq. " 
kinematic distance (Tableland 
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with our 



4.6) yields 



0.6 a% Gyr 



(19) 



at a radius of five arcminutes where d,2 is the distance 
to NGC 6397 divided by 2 kpc and we have depro- 
jected the velocity and proj ected stellar densities usin g 



an inverse Abel transform (Leonard & Mcrritt 



1989). 



If we use the velocities measured by MlViyi at five ar 
cminutes (the approximate radius of ou r field) with the 
stand ard candle distance of 2.53 kpc (Gratton et al. 
2003 ) yields the larger value of 



r r « 0.8 4. 5 3 Gyr. (20) 
We will revisit these timescales in the following section. 

4.5. Proper-Motion Distribution 

Our first focus will be the sample of about 2000 
main-sequence stars with 19.5 < F814W < 24.5 as 
in the previous section. As MAM06 did, we look 
at the differential distribution of proper motions in 
various directions as well as the total proper motion. 
These directional distributions are presented in Fig. [14} 
Each of the distributions has been fit with a Gaus- 
sian (Maxwcllian) distribution of proper motions in 
one dimension; the value of the standard deviation 
of the distribution is taken to be the value of a for 
the sample of proper moti ons. O ne can conclude from 
an Anderson-Darling test ( |1952 ) that the proper mo- 
tions in the radial and right ascension directions are 
not sampled from a Gaussian distribution at the 95% 
confidence level. The conclusion that the proper mo- 
tions in the tangential and declination directions are 
not drawn from a Gaussian exceeds 99% confidence. 
The values of a for the various directions do not differ 
significantly fro m e ach other as one would expect from 
the results of ; 



4.3 



We can directly compare the various distributions 
with each other and look for significant differences us- 
ing a KS test. Because the KS test is most sensitive 
to changes in the median values (here expected to be 
zero for all of the one-dimensional distribution), we 
compare the distribution of the absolute values of the 
proper motions. In particular the distribution of the 
proper motion toward or away from the cluster cen- 
ter (radial) or perpendicular to this direction (tangen- 
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tial) do not differ significantly as the KS test yieids a 
p— value of 57%. On the other hand, the hypothesis 
that the proper motions in right ascension and decli- 
nation are drawn from the same distribution can be 
rejected with nearly 99% confidence according to a KS 
test, indicating that there may be some residual sys- 
tematics in the proper-motion measurements; however, 
the values of the dispersions of the best-fitting Gaus- 
sian agree well within the ninety percent confidence 
regions, so the differences are more subtle than the 
width of the distributions. 



TABLE 4 

KS AND WlLCOXON PROBABILITIES 



T he t heoretical model that we developed in §<: 3.3 
and |4.1| exhibits an isotropic velocity distribution, so 
it makes most sense to compare the model with the dis- 
tribution of the total proper motion. We shall examine 
two theoretical models. First, the model yields an esti- 
mate of the proper-mot ion distribution over the entire 
field as depicted in Fig. |15| This p rope r-motion distri- 
bution results from convolving Eq. |A9| in the appendix 
with the window function in spherical radius in Fig. [3] 
This proper model distribution yields a KS probabil- 
ity of two percent for the null hypothesis that the stars 
from the prop er-m otion sample are drawn from it. The 
results of Fig. [13] indicate a second possible treatment. 
The figure depicts the model results for the stellar col- 
umn density and the dispersion (a). Although the fit 
to the proper motion appears good, the model does 
not track the observed column density well. In fact 
the observed column density does not actually decrease 
monotonically with radius, indicating that possibly the 
cluster is a bit clumpy. How ever, we can also convolve 
the distribution in Eq. |A10| with the number of stars 
as a function of radius in Fig. [13] This results in the 
"Scaled Model" in Fig. [15] that yields a KS probability 
of thirty-three percent for the null hypothesis that the 
observed proper motions of the nearly two thousand 
main-sequence stars in the proper-motion sample are 
drawn from the theoretical model. 

An examination of the differential distribution of the 



total proper motion yields additional details (Fig. 16 1. 
The width is similar to the one-dimensional results m 



Fig. 14 The distribution of stars cuts off much more 
quickly than the best-fitting Rayleigh distribution (the 
distribution of speeds if the velocity in each of two di- 
mensions follows a Gaussian) . The latter gives 49 stars 
with total proper motions exceeding one milliarcsecond 
per year; there were only six observed in the sample. 
If we assumed that the best fitting Rayleigh distribu- 
tion was the underlying true distribution the chance of 
finding six when we expect 49 is about 10~ 14 , assum- 
ing Poisson statistics. The cutoff at about one milliarc- 
second per year results because the cluster has a finite 
escape proper motion at this radius of ro ughly one mil- 
liarcsecond per year. The results of § 4.1 prod uce a dis- 
tribution of total proper motions using Eq. ( A9 ) con- 
volved with the window function shown as the "Model" 



curve in Fig. 16 ) or convolved with the observed num- 
ber of stars as a function of radius from Fig. |13| shown 
as the "Scaled Model" curve. 

To contrast the distribution of white dwarfs and 
main sequence stars, we examine the subsamples de- 
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Kolmogorov-Smirnov (above diagonal) and Wilcoxon (below di- 
agonal) probabilities in percent from comparing the proper- 
motion distribution of one of each pair of samples with the other. 
The BWDN values compare the BWD sample with the other 
samples with added Gaussian errors such that the root-mean- 
square error in other samples equals that of the bright-white- 
dwarf sample. 



lineated in Fig. [2] and Table [3] Fig. [17] depicts the cu- 
mulative proper-motion distributions for the four sub- 
samples, and Table[2]gives the KS and Wilcoxon proba- 
bilities for each pair. Lower probabilities indicate that 
the null hypothesis that the two samples are drawn 
from the same underlying distribution is unlikely. The 
stars in the bright WD sample and faint MS sample 
span the same magnitude range, but a comparision 
of their proper-motion distributions yields p— value of 
0.02 by the KS test. The main sequence stars with 
F814W « 19 (middle MS) are likely to have similar 
masses to those of the brightest white dwarfs. These 
stars yield a larger p— value of about 0.12. 

The bright main-sequence stars yield the largest 
p— value in a comparison with the bright white dwarfs 
of 0.22. The null hypothesis that the bright white 
dwarfs and the bright main-sequence stars are drawn 
from the same distribution is the most difficult to re- 
ject by these tests. If we suppose that these two groups 
of stars indeed follow the same underlying proper- 
motion distribution, we would conclude that the bright 
WD stars have not yet relaxed to reflect their smaller 
masses. The observed luminosities of these stars pro- 
vide constraints on their ages. The median magnitude 
of the bright-white-dwarf s ample is 24.27 in I— band 
with the transformations of Sirianni et al.| (|2005|) from 
F814W and F606W to I. Using the w Eite^dwarf 
distance of 2.53 kpc (Hansen ct al. 2007), a white 



dwarf mass of 0.53 M ( |Moehler et ah 2004 K alirai 
let al.| | 2009[ ) and the cooling models ot Bergeron et al. 



(|1995p yields a median age of the bright WD sa mple 



of 0.5 Gyr. The shorter kinematic distance (see § 4.6 1 



yields a median age of 0.9 Gyr. We will simply take 
the mean of these two determinations yielding a me- 
dian age of about 0.7 Gyr. 

If the white dwarfs with a median age about 0.7 Gyr 
have a similar proper-motion distribution to their pro- 
genitors even though their masses are substantially 
less, these stars have not yet relaxed to energy equipar- 
tition, and this fact yields a direct estimate of the relax- 
ation time of the cluster in the ACS field. This conclu- 
sion agrees with dynamical estimate of the relaxation 
time for the white dwarfs from Eq. [19] of 1 Gyr at the 
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Fig. 14. — Differential distributions of the observed proper motions in various directions and best-fitting Gaussians. The values of 
a are given in milliarcseconds per year. 



median projected radi us of the sub-sample, five arcmin- 
utes (Fig. pi). Fig. 12 apparently shows that the older 
white dwarfs have a similar velocity dispersion to the 
main-sequence stars. The median F814W magnitudes 
of the next two fainter white-dwarf points are 25.66 
and 26.35, corresponding to white-dwarf ages of abo ut 
2.5 and 5.2 Gyr respectively (iBergeron et al. 19951, a 



factor of three to seven older than the bright WD sam- 
ple. This yields a loose up per limit to the relaxation 
time in agreement with Eq. |19[ however, the errors in 
proper motions also begin to dominate at these mag- 
nitudes, making it difficult to draw conclusions as we 
must rely more strongly on the subtraction of the ob- 



servational errors. Better constraints on the proper 
motions at faint magnitudes from future observations 
could allow us to place a more reliable upper limit on 
the relaxation time. 

The proper-motion errors increase at faint magni- 
tudes because the poisson noise in the pixels used to 
determine the position increases as the flux decreases. 
Brighter than F814W = 19.5 saturation begins to set 
in for the deep images, so a smaller set of shorter ex- 
posures is used for the astrometry of brighter stars. 
The proper-motion errors increase abruptly here and 
again decrease with brighter magnitudes. This results 
in a variation of the typical proper motion error for 
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Fig. 15. — Cumulative distribution of total proper motion for 
the main-sequence stars, the basic model and a model scaled to 
yield the observed column density in Fig. |13| 
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Fig. 16. — Differential distributions of total proper motion for 
the main-sequence stars, the basic model and a model scaled to 
yield the observed column density in Fig. |13| and the best-fitting 
Gaussian. The value of a is given in miliiarcseconds per year. 

the various samples as shown in Table [3j These er- 
rors contribute to the observed dispersions so they 
hav e be en subtracted from the dispersions in Fig. [12] 
and |13| They also b road en the cumulative distribu- 
tions depicted in Fig. |17| and bias the probabilities in 
Table [4] From Table |3j it is apparent that the typical 
error in the white dwarfs is substantially larger than 
for the other samples and that the errors in the other 
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Fig. 17. — Cumulative distributions of total proper motion for 
the subsamples. The observed distributions are given by the 
solid curves, and the model distributions using the masses in 
Table [3] are traced by the dashed curves 

samples are nearly equal to each other; therefore, we 
artificially increase the proper-motion errors in these 
samples by adding normally distributed random num- 
bers to each of the proper-motion measurements, so 
that all four distributions have the same root-mean- 
square error as the bright white dwarfs. Although 
the faint-main-sequence sample (FMS) and the white- 
dwarf sample cover the same range of apparent magni- 
tude, the white dwarfs are typically fainter within that 
range than the main-sequence stars (see Tab. [3]) and 
hence exhibit larger proper-motion errors. 

We then calculate the KS and Wilcoxon probabil- 
ities for the new subsamples against the white-dwarf 
distribution. One thousand new subsamples are gen- 
erated in this way, and the average probabilities are 
tabulated. These results are depicted in the column 
and row labelled "WDN" in Table [4] By making the 
distributions even broader than the white dwarf dis- 
tribution, the additional error decreases the p— values 
that result from comparing one sample to another, ver- 
ifying that the differences in the distributions do not 
result from differences in the errors. 

The distributions depicted in Fig. [17] merit further 
comment. The faint MS sample appears to have typ- 
ically larger proper motions than the more massive 
middle MS sample, and the brightest MS stars have 
slightly smaller values still. This is in agreement with 
the values of a in Table [3] and with the expectations 
of equipartition. However, from Fig. [5j [12] and [17] we 
can seen that the velocity dispersions within our field 
expected from the models vary little with mass (in fact 
less than is observed), but observations nearer to the 
core would be a better probe of the dynamical state of 
the various stars in the cluster. 



18 



4.6. The Kinematic Distance of NGC 6397 

Fig. [5] depicts the observed proper- motion measure- 
ments as a function of distance from the center of 
the cluster along with the measurements of the radial- 
velocity dispersion converted to proper motions at a 
fiducial distance of 2.55 kpc ( Richer et al. 2008) and 
our best-fi tting distance of 2 kpc. |Leonard fc Mer- 



ritt 



( 19891) argue that for a spherical stellar cluster the 



proper motions themselves will reveal t he ha llmark of 
an anisotropic velocity ellipsoid. In § 4.3 we argue 
that there is no compelling evidence for anisotropy ei- 
ther in the sample as a whole or in various subsamples; 
therefore, it is natural to compare the observed proper- 
motion dispersions with the radial velocity dispersions 
to obtain a kinematic (or geometric) estimate of the 
distance to the cluster. 

The best-fitting distance with the smallest statistical 
error in Table [5] is obtained by estimating the proper- 
motion dispersion of all the main-sequence stars in the 
proper-motion sample, 0.369 ±0.009 mas/yr. We then 
fit the MM91 data linearly with projected radial dis- 
tance to estimate the radial-velocity dispersion at the 
median distance of the stars in our sample (5.16 ar- 
cminutes) yielding 3.12 ±0.11 km/s. This yields a dis- 
tance of 1.78 ± 0.07 kpc. Because it involves the most 
data, this estimate has the smallest statistical confi- 
dence interval. However, there are two possible sources 
of systematic error that we would like to address. First 
MM91 actually measured the radial- velocity disper- 
sion at 5.22 arcminutcs near the median radius of our 
dataset to be 3.3 ± 0.8 km/s. Combining this velocity 
with the dispersion of the entire dataset yields a larger 
distance of 1.9 ± 0.5 kpc. Second, the stars studied by 
MM91 were among the brightest in the cluster so it is 
natural to compare the brightest stars in our sample 
with their data. Fig. [12] shows that the proper-motion 
dispersion decreases slightly with increasing flux. The 
dispersion (a) of the brightest 200 stars in our sam- 
ple is 0.33^q'q3 mas/yr with a slightly smaller pro- 
jected median distance of 5 arcminutes. To obtain 
a more robust kinematic estimate of the distance to 
NGC 6397 we have used the unpublished radial veloci- 
ties from MM91 (Meylan, private communication). We 
restrict the MM91 dataset to 46 stars with projected 
distances from the center of NGC 6397 between 3 and 
7.2 arcminutes; the median projected distance of this 
sample from the center is slightly smaller at 4.5 ar- 
cminutes. The dispersion of this bright overlap sample 
is 3.5^q g km/s after correcting for measurement un- 
certainties, yielding a distance estimate of 2.2^7 kpc 
with a ninety-percent confidence interval. We obtain 
our best distance estimate by fitting the MM91 data 
in radius and comparing the dispersion to that of the 
brightest 200 stars in our sample (1.9 ± 0.2 kpc). Ad- 
ditionally we correct for the difference in the median 
mass of these two samples (O.84M for the giants and 
O.74M for the bright proper motion sample). This 
yields a distance estimate of 2.0 ± 0.2 kpc and bal- 
ances biases against statistical errors. For reference all 
of these distances are presented in Table [5] 



TABLE 5 
Distance Estimates to NGC 6397 



Method Value [kpc] 

Standard-Candle Estimates 
Subdwarf fit in B and V 2.53 ± 0.05 



References 



Subdwarf fit in b and y 


2.58 


± 


0.07 








Subdwarf fit 


2.67 


± 


0.25 


Reid & 


Gizis 


1998 


White-dwarf fit 


2.55 


± 


0.11 


Hansen 


et ai. 


2007 



Kinematic Estimates 



MS sample (radial fit) 
MS sample (at 5.22') 
Bright sample (fit) 
Bright sample (fit/corr) 
Bright overlap sample 



1.78 ±0.07 
1.9 ±0.5 
1.9 ±0.2 
2.0 ±0.2 
2.2 4 



MM91/Here 



,+0.5 
-0.7 



These successive distance estimates trade larger sta- 
tistical errors for reduced systematic errors that arise 
from the heterogeneity of the stars and through in- 
terpolating the velocity profiles in radius. The final 
estimate of 2.2^q7 kpc compares stars of nearly the 
same mass at the same median radius, so it minimizes 
the systematic errors; however, it uses only a portion 
of both samples, so the statistical errors are necessar- 
ily larger. Increasing the size of the sample of stars 
with radial-velocity measurements could dramatically 
decrease the uncertainities of this kinematic distance 
estimate. 

4.7. The Mass of NGC 6397 



Leonard & Merritt (1989) outlined several mass es- 
timators tor the open cluster M35 using the observed 
proper-motions. Among these the most straightfor- 
ward is 



(GM r ) = 



16 



R 



1 



(21) 



where vr is the component of the proper-motion to- 
ward the center of the cluster and vt is the tangential 
component of the proper-motion. This estimator yields 



8.3 ±0.4 x 1O 4 4 53 M 



(22) 



where we have summed over all of the stars along the 
cluster main sequence and white-dwarf cooling tracks 
and excluded those stars with proper motion errors 
greater than 0.4 mas/yr or proper- motions greater 
than 5 mas/yr. The confidence interval gives a ninety- 
percent confidence region obtained through bootstrap- 
ping the sample. 

We use the power-law fits to the column density and 
velocity dispersion (Eq. 16 ) of the cluster to deproject 
the proper-motion dispersion, density and their gradi- 
ents to obtain an estimate of the mass enclosed within 
the outer radius of our field. We can simply apply 
Jeans equation 



d In ai 



din? 



dlnr dlnr 



GAL 



(23) 



to obtain an estimate of the enclosed mass within a 
given spherical radius (M r ). The symbols o~ r and n 
denote the radial velocity dispersion and the number 
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density of stars. In this equation v r and j3 quantify 
the rotation of the cluster and the anisotropy of the 
velocity distribution; we see evidence for neither in our 
data, so we assume these vanish to yield 



M r 



8.7 ±0.6 x 10 4 d: 



4j3 



2.53 



0.6 



M. 



(24) 



Whether one uses a or a in Eq. 
mass estimate (Eq. 24 1 within t 



23] does not affect the 
We used 



;he errorbars 

proper motions and positions themselves without bin- 
ning in radius to determine the best fitting power-law 
relations (Eq. 171 for the column density and velocity 
dispersion, yielding a mass estimate of 



M r = 8.5 ± 0.6 x 10 4 <^ 53 



0.64 



M 



0- 



(25) 



We can combine these estimates with the light within 
this projected radius to yield 

L V<R = 4.0 x 10 4 4 53 L y , o (26) 

a nd t he model for the light distribution outlined in 
§ |4.1| allows us to estimate the luminosity within the 
spherical radius at 



L 



V,r 



3.5 x 10 4 4 53 L y , o 



(27) 



Scaling the mass within seven arcminutes (Eq. 25 1 by 
the ratio of the total ligh t (Table [TJ to that within 
seven arcminutes (Eq. [27] ) yields an estimate of the 
total mass of the cluster of 



M = 1.1 ±0.1 x 10 5 £ „M 



2.53 1V1 



(28) 



in a greem ent with the model mass (Eq. 12 ) found 
in § |4.1| This gives a mass-to-light ratio of 2.4 ± 
O.Sc^Mq/L^©; including the uncertainty in the dis- 
tance yields 2.4 ± O.SMq/Lv,©- This agrees with the 
previous result of MM91 of 2.1 ± 0.1 at an assumed 
distance of 2.4 kpc. The mass-to-light ratio of the to- 
tality of the stellar population of NGC 6397 range from 
1. 4 to 2.0 accordi ng to MM91 and about 1.2 according 



to Drukier ( 1995[ ) . Our present mass-to-light measure- 
ment does not require any dark contribution to the 
mass of NGC 6397. 

4.8. Stellar Escapers 

Fig.[5]also depicts the escape velocity as a function of 
distance from the center of the cluster. The escape ve- 
locity (here given as an escape proper-motion) is about 
three times the velocity dispersion for the range of radii 
probed by the observations. This is the escape veloc- 
ity at a given distance (not projected) from the center 
of the cluster; therefore, a star with a proper-motion 
larger than this at a particular projected radius will 
escape. Exceeding this proper-motion is a sufficient 
condition for escape. Other stars whose actual three- 
dimensional distances are further from the cluster cen- 
ter may also escape as can stars with significant motion 
along the line of sight. For an isotropic velocity distri- 
bution the three-dimensional velocity is typically only 
about twenty percent larger than the proper-motion; 




Ha [mas/yr] 



Fig. 18. — Proper motions of the objects identified as stars 
in the ACS field and lying near the main sequence and white 
dwarf cooling track of NGC 6397. The mean proper-motion of 
the cluster members has been subtracted. The red circle denotes 
twice the escape proper-motion from the inner edge of the ACS 
field. 



therefore, we will use a conservative criterion and con- 
sider stars that exceed the escape proper-motion as 
potential escapers. 

Fig. [18] depicts the proper motions of all stars along 
the main-sequence and white-dwarf tracks. The tight 
cloud centered near zero proper motion are the clus- 
ter members whereas the more diffuse cloud at /i Q rs 
— 12 mas/yr consists of field stars that happen to have 
fluxes similar to cluster main-sequence or white-dwarf 
stars. A comparison with the proper motions of all 
stellar sources in Fig. [l] demonstrates that the flux cri- 
teria removes nearly all the field stars while affecting 
the potential cluster stars more modestly. We define 
the region of potential escapers as having a total proper 
motion between once and twice the escape proper mo- 
tion for stars at the inner edge of the ACS field. This 
upper limit is indicated by the red circle. Looking at 
the density of presumably field stars outside this circle, 
one expects fewer than one field star within the narrow 
range of color and proper motion. 

Table [6] lists those stars that lie on or near the clus- 
ter main-sequence or white-dwarf tracks whose total 
proper motion exceeds the escape proper motion at 
the tangent point for the projected radius of the star. 
We have verified that all of these potentia l escapers 
satisfy the im age-quality criteria outlined in [Anderson 
et al. (2008). Fig. 19 depicts the ratios of the stellar 
proper motions to the escape proper motion (Tab [6] 
lists those stars between the two circles). Looking at 
the proper motions themselves in Fig. 19 the potential 
white-dwarf escapers (green) appear to be a smooth 
extension of the rest of the stars; whereas the poten- 
tial main-sequence (red) escapers appear to lie beyond 
the tail of the main-sequence distribution. The posi- 
tions of the potential escapers in Fig. [20] do not appear 
extraordinary; however, th eir locations on the color- 
magnitude diagram in Fig. 21 are. The potential es- 
caping white dwarfs are among the faintest stars with 
measured proper motions in the sample and therefore 
are prone to greater error. This agrees well with the 
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TABLE 6 

Cluster stars with proper-motions that exceed the local escape proper motion 



ID 


F814W 


F606W-F814W 


Mtot 
[mas 


,fJ>e 

yr- 1 ] 


C12000 

[arcsec] 


[arcsec] 


[mas vr 1 

L J I 


[mas yr" 


j 


Up 

[mas vr 1 

L J J 


Comment 


14960 


18.92 


1.03 


1 


27. 


1 


10 


40.38 


84.95 


1.24+0.08 


-0.28+0 


08 


0.77 ± 0.08 


MS, short 


35641 


18.93 


1.00 


1 


52. 


1 


15 


33.10 


167.75 


1.16+0.05 


-0.98+0 


05 


1.00 ± 0.05 


MS, short 


39398 


18.93 


1.02 


1 


32. 


1 


13 


71.20 


180.95 


1.17±0.08 


-0.62+0 


08 


1.00 ± 0.08 


MS, short 


34104 


19.04 


1.04 


1 


11. 


1 


15 


33.82 


162.30 


1.30+0.05 


-0.56+0 


05 


0.95 ± 0.05 


MS, short 


16272 


19.27 


1.08 


1 


12. 


1 


10 


35.40 


91.00 


0.81+0.03 


0.76+0 


03 


-0.11 ± 0.03 


MS, short 


25395 


21.12 


1.38 


1 


34 


1 


11 


70.85 


129.75 


1.26+0.01 


0.44+0 


01 


0.48 ± 0.01 


MS 


14687 


23.42 


1.98 


1 


11. 


1 


06 


134.25 


83.70 


-0.32+0.05 


-1.06+0 


05 


0.41 ± 0.05 


MS 


19308 


25.67 


0.74 


1 


10. 


1 


08 


96.95 


102.95 


-0.16+0.21 


-1.09+0 


21 


0.56 ± 0.21 


WD 


21043 


26.14 


0.91 


1 


32. 


1 


10 


68.60 


112.45 


-0.74+0.35 


1.09+0 


35 


—0.98 ± 0.35 


WD 


27645 


26.18 


0.96 


1 


19, 


1 


07 


137.35 


139.00 


-1.18+0.27 


-0.14±0 


27 


-0.77 ± 0.27 


WD 


34340 


26.23 


0.93 


x 


26. 


1 


11 


90.70 


163.20 


-1.10+0.33 


-0.61±0 


33 


-0.46 ± 0.33 


WD 


9666 


26.29 


0.96 


1 


21. 


1 


08 


48.60 


59.70 


1.07+0.32 


-0.56+0 


32 


0.89 ± 0.32 


WD 


15213 


26.30 


1.01 


1 


12. 


1 


01 


184.80 


86.15 


-0 95+0 36 


0.59+0 


36 


— 1 00 + 36 


WD 


4333 


26.39 


1.07 


1 


29. 


1 


08 


10.77 


30.96 


-0.93±0.33 


-0.90+0 


33 


0.30 ± 0.33 


WD 


28585 


26.40 


1.05 


1 


32. 


1 


11 


23.43 


142.55 


1.14+0.36 


-0.66+0 


36 


0.95 + 0.36 


WD 


25008 


26.44 


1.03 


1 


45, 


1 


05 


181.60 


128.10 


0.12+0.26 


-1.44+0 


26 


0.57 + 0.26 


WD 


3445 


26.54 


1.14 


1 


37, 


1 


07 


20.72 


25.39 


0.55±0.31 


1.25+0 


31 


-0.62 + 0.31 


WD 


30992 


26.55 


1.12 


1 


30. 


1 


07 


141.70 


151.00 


0.94±0.28 


-0.90+0 


28 


0.97 + 0.28 


WD 


433 


26.61 


1.22 


1 


35, 


1 


03 


123.15 


6.78 


0.63±0.37 


1.19+0 


37 


-0.38 + 0.37 


WD 


27465 


26.62 


1.18 


1 


55, 


1 


08 


120.70 


138.20 


1.19±0.35 


1.00+0 


35 


0.27 + 0.35 


WD 


31935 


26.65 


1.16 


1 


30, 


1 


15 


24.08 


154.50 


0.81+0.28 


-1.02+0 


28 


0.99 + 0.28 


WD 


16969 


26.67 


1.20 


1 


52. 


1 


10 


42.69 


93.90 


-0.57±0.34 


-1.40+0 


34 


0.49 ± 0.34 


WD 


30725 


26.72 


1.27 


1 


42, 


1 


06 


176.30 


150.00 


0.54±0.30 


1.31+0 


30 


-0.09 + 0.30 


WD 


25724 


27.00 


1.26 


1 


65, 


1 


06 


163.30 


131.10 


0.26+0.39 


1.62+0 


39 


-0.38 + 0.39 


WD 



The quantities c*2000 an d (52000 gi ve the offset in arcseconds from the origin of the coordiates of the field at a = 17 h 40 m 56!72 and 
5 = -53°45'36'!8. 

apparent smoothness of their proper-motion distribu- 
tion with the rest of the white-dwarf population. 

On the other hand, the location of the potential 
main-sequence escapers in Fig. |2l] presents a bit more 
of a puzzle. All but two of these stars lie between 
18.9 < F814W < 19.3 where the astrometry begins to 
use the shorter exposures in which these stars are not 
saturated. In fact in these short exposures the proper- 
motion errors are comparable to those of F814W « 24 
in the longer exposures and approach 0.1 mas/yr. Al- 
though these errors are not as large as those for the 
potential white-dwarf escapers, there are many more 
stars at these magnitudes, so one would expect a num- 
ber of outliers. The final potential escapers lie at 
F814W 21 and 23 where the proper- motion errors 
are typically less than 0.05 mas/yr. Statistically it 
is unlikely that these objects are field stars for the 
reasons mentioned earlier. The brighter, faster mov- 
ing star (ID 25395) with a total proper motion of 
1.34 mas/yr appears cle arly beyond the tail in the dis- 
tribution shown in Fig. [16] From Tab. [6] we see that 
its total proper motion exceeds its escape proper mo- 
tion by sixteen times its proper-motion error. On the 
other hand the fainter, slower moving star (ID 14687) 
exceeds its escape proper motion by about the value 
of its proper motion error, so its assignment as an 
escaper is less certain. In fact given the number of 
stars slightly below the escape threshold it is likely 
that this star is not an escaper but rather measure- 
ment errors pushed it above the threshold; therefore, 




Fig. 19. — The proper motions of potential stellar escapers. 
Only objects with proper-motion errors less than 0.4 mas/yr are 
included (the minimum escape proper motion is 1.3 mas/yr). 
There is no proper-motion cutoff. The arrow indicates the di- 
rection of the center of the cluster. The proper motions of stars 
along the main sequence and white-dwarf track are depicted with 
red triangles and green squares respectively. The most likely es- 
capers are circled. The escape proper motion is denoted by the 
solid red circle and twice its value by the dashed blue circle. 
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Fig. 20. — The positions of potential stellar escapers. The 
center of the cluster is located beyond the upper-left corner of 
the field. The red dots show the positions of all stars in the 
sample, the potential escaping MS stars are blue triangles, and 
potential escaping WD stars are green squares. The motions of 
two most likely escapers over the next fifty thousand years are 
depicted. 




F606W-F814W 

Fig. 21. — The color and magnitude of potential stellar esca- 
pers. The red dots show the positions of all stars in the sample, 
the potential escaping MS stars are blue triangles, and WD stars 
are green squares. 



we estimate that there is one escaper out of a sample 
of 3,245 stars. Typically a star on an escape trajectory 
will leave the cluster within a crossing time, r c ~ r/v « 
5'/0.3 mas yr" 1 « 10 6 yr. This yields an estimate of 



the evaporation timescale, (dlnN/dt) -1 ~ 3 Gyr, or 
about three relaxation timescales. This evaporation 
timescale does not indicate that the cluster will van- 
ish over the next 2 Gyr. On the contrary the escaping 
stars as apparent from Fig. 



21 arc less massive than 



average so the mass-loss timescale will be perhaps two 
or three times longer; furthermore, this timescale is 
the time over which the cluster loses about two-thirds 
of its stars, not all of them. This agrees with the de- 
struc tion timescale calculated by (Gnedin & Ostriker 
1997) of about 4 Gyr. This star is moving away from 
the centre of the cluster, although it does not travel 
on a strictly radial trajectory. This is not surprising 
as the diffusion of stars in velocity to become unbound 
is not dominated by a few strong encounters but by 
many small ones. The final of these small encounters 
could have occurred anywhere within the cluster not 
necessarily near the centre, so the trajectory of a star 
as it escapes need not be radial. 

5. DISCUSSION 
5.1. Dynamics 

The proper motion measurements by ACS in 
NGC 6397 reveal gross agreement with the features 
of a mass-segregated lowered isothermal distribution 
function for the stars. On the face of it this conclu- 
sion is not particularly surprising because such a model 
was conceived to describe globular clusters such as 
NGC 6397. There is no strong evidence for anisotropy 
in the proper motion distribution, neither in the stellar 
population as a whole nor in subsamples of stars of par- 
ticular apparent magnitude or at particular projected 
distances from the cluster center. 

Three discrepancies with the lowered isothermal dis- 
tribution are apparent. First, the effects of mass seg- 
regation appear larger in the data than predicted in 
the model under the assumption that m,^ 2 is con- 
stant. The velocity dispersion in the ACS field de- 
pends modestly on the apparent magnitude (the mass) 
of the stars, as expected for the model; however, the 
mass dependence appears to be slightly stronger in the 
data than in our model (Fig. 12). Furthermore, the ra- 



dial distribution of stars in our sample depends more 
strongly on mass than in the models (Fig. [6| ; these dif- 
ferences are stronger than we expect from incomplete- 
ness (Fig. p}. Second, the proper motion distribution 
of young white dwarfs (younger than 0.7 Gyr) appear 
to not satisfy equipartition with the other stars of the 
cluster. Their proper motions are probably smaller 
than those of stars of similar masses, and their proper 
motions are most similar to the motions of the most 
massive stars observed in the cluster, i.e. stars similar 
to their progenitors. The young white dwarfs do not 
appear to have relaxed, yielding the first direct mea- 
surement of the relaxation time in a globular cluster of 
about or greater than 0.5 Gyr. This value agrees with 
the theoretical expectations from the number of stars 
in NGC 6397 near the half-light radius and the velocity 
dispersion of the cluster. Thirs, there are some stars 
with sufficient energy to escape the cluster whereas 
the lowered isothermal distribution function contains 
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no escaping stars. The abundance of these stars along 
with the observed crossing time yields an estimate of 
the evaporation time of the cluster of about 3 Gyr. 
Because the typical escaping stars have masses smaller 
than average, the mass-loss timescale is longer perhaps 
by a factor of a few. This constraint agrees with the- 
oretical expectations (Gnedin & Ostriker 1997) of a 
destruction timescale of about 4 Gyr. 

5.2. Properties 

Because the ACS field lies beyond the half-light ra- 
dius of the cluster, the measurements yield an esti- 
mate for the total mass of the globular cluster with 
only modest extrapolation. The resulting mass is 
1.1 ± 0.1 x lO 5 ^ 53M0 in agreement with the result 
of MM91. Proper motions, unlike radial velocity mea- 
surements, probe the isotropy of the velocity distri- 
bution, so they give a m ore robust mass estimate 
(Leonard & Merritt 1989). Because the proper mo- 
tions show no evidence for anisotropy, the line-of-sight 
velocity measurements of MM91 who assumed isotropy 
gave a similar mass estimate and mass-to-light esti- 
mate. This mass-to- light ratio (of about 1.5 in solar 
units) does not indicate a need for dark matter in the 
cluster to explain the observed dynamics. 

The major difference with previous results is the 
kinematic estimate of the distance to NGC 6397. 
These estimates range from 1.78 ± 0.07 kpc to 
7 kpc as summarized in Table [EJ The larger 
estimates mitigate against various known biases at 
the expense of increased statistical uncertainties. Our 
best estimate with the current data of 2.0 ± 0.2 kpc 
with ninety-percent confidence was obtained by fitting 
a linear model to the radial velocity dispersion mea- 
sured by MM91 and comparing this with the disper- 
sion of the proper motions of the brightest stars in our 
sample and by correcting for the difference in mass 
between the bright-star sample and the giants. Al- 
though our longest distance (and most uncertain) es- 
timate does agree with the previous standard-candle 
estimates, the most precise of which is 2.53 ± 0.05 kpc 
within their individual ninety-percent confidence inter- 
vals, there is some tension between our best estimate 
(2.0 ± 0.2 kpc) and the standard-candle estimates. A 
short kinematic distance to NGC 6397 is similar to the 
short distance found by MAM06 for 47 Tuc that was 
also about twenty percent smaller than the standard 
candle estimates. This kinematic distance estimate is 
still plagued by potential systematic errors because the 
stars with proper motions are not the stars with radial 
velocities; therefore, some extrapolation over position 
in the cluster and stellar mass is required to get the dis- 
tance estimate. Such a short distance, if indeed cor- 
rect, could pose a challenge for stellar evolution and 
white dwarf cooling models; however, it does indicate 
a need to obtain more precise and model-independent 
measurements of distances to globular clusters. 

5.3. Future Directions 

In light of the data presented, further modelling 
of the cluster is warranted. In particular a larger 



N— body model than that recently performed by |Hur-| 
ley et al. ( 2008 1 to account for both the observed star 



count and the observed proper motions could dramat- 
ically increase our understanding of this well-studied 
cluster. In particular the improved model of the clus- 
ter should include the effects of disk sh ocking to con- 
sider the effects of the observed orbit (Kalirai et al. 
2007) of the cluster. As mentioned in the introduc- 
tion, the number of the stars in the outskirts of the 
Hurley et al. (20081 model were insufficient for a de- 
taiied comparison with the data. A model resulting in 
a final cluster with more than 150,000 stars would be 
sufficient for a comparison in this field. An observa- 
tional study closer to the core of NGC 6397 would be 
easier to compare with numerical models because the 
higher stellar densities both in the models and the data 
may allow a detailed comparison with current models. 
Combination of the outer and the inner observations 
for this cluster or 47 Tuc would be important. Our 
attempts to obtain a kinematic estimate of the dis- 
tance to NGC 6397 have been thwarted by potential 
systematic errors in particular, a lack of overlap be- 
tween the stars with proper motions (presented here) 
and those with radial velocities (MM91). This obsta- 
cle also affected the work of MAM06 for 47 Tuc. A 
natural way forward would be to obtain radial velocity 
measurements for the brightest stars in our proper mo- 
tion sample with F814W ~ 16. This would yield the 
first significant sample of stars with known velocities 
in all three directions for a globular cluster and would 
prove a powerful tool in probing the cluster as well as 
obtaining a distance estimate. 

Given the relative vicinity of NGC 6397 to Earth 
(and 47 Tuc as well for that matter) and the high- 
precision astrometry of ACS, a natural next step would 
be to measure the distances to these clusters with par- 
allax. Both astrometric epochs for this work were ob- 
tained at the same time of year, removing the effects 
of parallax to fo cus on the proper m otion of the clus- 
ter and its stars. Kalirai et al. (2007) found a sufficient 
number of background galaxies to fix the astrometry to 
a precision of about 0.2 mas. Even with ACS observa- 
tions at a different time of year, this level of precision 
is insufficient to improve on the current precision of 
the measurement of the distance to NGC 6397. Two 
potential ways forward would be to find another way 
to fix the astrometry or to obtain proper motions and 
radial velocities for the same stars. In the more distant 
future GAIA should be able to measure the distance 
to NGC 6397 and other nearby clusters. Distances to 
these nearby metal-rich and metal-poor clusters would 
have dramatic implications beyond the study of glob- 
ular clusters themselves from detailed stellar models 
to the extragalactic distance ladder, extending to the 
edge of the observable Universe. 
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APPENDIX 

The observations presented here probe beyond the half-light radius of NGC 6397. In particular a useful approx- 
imation is the limit where the value of v e <C o as is appropriate in this region. In this limit we have 



J n 



1 



„"+2 



,,n+5 



dv 



1 Ve 

a 2 (n + 3)(n + 5) 



(Al) 



so ratios such as (v 2 ) — J2/J0 no longer depend on a. This particular limit applies near the edge of the cluster 
regardless of the underlying distribution function as long as the distribution function is approximately linear in 
energy near the escape threshold. Furthermore, the linear approximation of the distribution function underesti- 
mates the number of stars with small velocities; therefore, it will overestimate the value of higher moments of the 
distribution function relative to lower ones. In particular 



(v 2 ) = y £ 



J 1DI 
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(A2) 



yielding an estimate of the escape velocity in terms of the local one-dimensional velocity dispersion. Let us define 
the line-of-sight integrated quantities 
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sec 2 9d9. 



Let us assume that the bulk of the matter lies within the radius R, so 

»- 2 GM GM GM 
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where cosf 



R/r t , so in this limit 
K n (R) « 
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If R <C r t we get 
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(A6) 



v e (R) « I ^M^42(v 2 D ) LOS * 3.1^1 



_D/LOS- 



(A7) 



As the ratio of the projected radius to the tidal radius increases, v e (R) approaches \/8y/ (w 2 i) )los- This not 
only provides a useful empirical estimator of the local escape velocity to find stars that may be escaping but 
also demonstrates that the observed velocity dispersion integrated along the line of sight for the outer regions 
of the cluster is independent of the value of cr; therefore, it is also ind epen dent of the mass of the stars even if 



Al 



and 



A6 shows that the density of stars in the 



equipartition holds and a cx to -1 / 2 . Further examination of Eq. 

outskirts of the cluster is proportional to a~ 2 but the radial dependence is independent of the value of tr; therefore, 
within a field where v e <C o we expect that the signatures of mass segregation will be weaker than regions where 

v e > cr- 

This paper looks at the distribution of the magnitude of proper motions that results from integrating the 
phase-space density over the allowed range of velocities along the line of sight, 
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where p\ is a normalizing factor. Of particular interest is the cumulative distribution of the magnitude of the 
proper motion of the stars. The number density of stars with proper motions greater than a given value v is 



n(> v\i 



Pi 



exp 




(A9) 
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One should note that one can get Eq. A9 simply by replacin g v e in Eq 
therefore, in the limit as the minimum velocity v vanishes, Eq. |A9| yields Eq 
course the velocity distribution of the sample along the line of sight 



N{> v\R) 



n{> v\r) 



rdr 



R 



R? 



of the main text) with y«| — v 2 ; 
as expected. The observable is of 



(A10) 



which must be calc ulate d numerically. Alternatively we can integrate n(> v\r) against the spherical window 
function defined in § 3.2 to obtain the proper-motion distribution for the entire field. 



